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ABSTRACT 

We compare calculated intensities of lines of C n, N I, N n O I and n with a published 
deep spectroscopic survey of IC 418. Our calculations use a self-consistent nebular 
model and a synthetic spectrum of the central star atmosphere to take into account 
line excitation by continuum fluorescence and electron recombination. We found that 
the Nil spectrum of s, p and most d states is excited by fluorescence due to the low 
excitation conditions of the nebula. Many C II and O II lines have significant amounts 
of excitation by fluorescence. Recombination excites all the lines from f and g states 
and most O n lines. In the neutral-ionized boundary the N I quartet and O I triplet 
dipole allowed lines are excited by fluorescence, while the quintet O I lines are excited 
by recombination. Electron excitation produces the forbidden optical lines of O I and 
continuum fluorescence enhances the N I forbidden line intensities. Lines excited by 
fluorescence of light below the Lyman limit thus suggest a new diagnostic to explore 
the inner boundary of the photodissociation region of the nebula. 

Key words: atomic processes- planetary nebulae: individual: IC 418 



1 INTRODUCTION 

With an apparent simple geometry the low excitation 
planetary nebula IC 418 lends itself to simpler models 
with fewer assumptions on its structure and kinematics 
than are generally necessary to reproduce observations 
of planetary nebulae (PN e). A recent deep sp ectroscopic 
survey of this object by ISharpee et alJ (|2003h (hereafter 
SWBH) identified several hundred li ne s in IC 418 from 
which ISharpee. Baldwin fc Williams! (|2004h (hereafter 
SBW) derived ionic abundances for CNONe ions rel- 
ative to H. Many dip ore-allowed line intensities from 
those ions in that survey cannot be explained by the 
recombination theory with the same abundances derived 
from collisionally excited lines (CELs). The intensities 
of recombination lines are often used to measure ionic 
abundances because they depend more weakly on tem- 
perature and density than CELs. Abundances derived 
from recombination lines, however, are usually higher 
than those derived from CELs in the Orion nebula 
llEsteban et all 1 19981: iBaldwin et alJ l2000l: lEsteban et all 
12004 iMesa-Delgado. Esteban. fc Garria-Roias! l200Sf ). 



many H 11 regions dc arcra Roias fc Esteban! 120071) 
and planetary nebulae (|Liu et al.l 19951. 12001b 



Luo. Liu fc Barlowll200ll:lLiu et al.ll2006l:lTsamis et al.ll2008l 
Garcia-Roias. Pena M. fc Peimbertl 120091 : lOtsuka et al.l 
20091 ). IWilliams et alJ i|200Sl ) showed that CELs probably 
give the correct abundances because the values derived 
from CELs are consistent with those derived from UV 
absorption lines in a sample of four PNe. The discrepancies 
of less than a factor of 2 found by SBW in IC 418 between 
abundances from CELs and dipole-allowed lines are smaller 
than in other objects. Nevertheless it is important to 
ensure that the intensities of the dipole-allowed lines are 
not enhanced by additional excitation mechanisms besides 
recombination whene ver they are used in tem perature or 
density diagnostics (|Fang. Storey fc Liu! |201ll ) (hereafter 
FSL) or abundance determinations. 

The usual procedure to determine a species X abun- 
dance directly from optical line intensities is to measure the 
ratio of the intensity of a line emitted by the species with 
respect to the H/3 intensity and assume that the ratio of col- 
umn densities producing the lines is equal to the abundance. 
Specifically, the abundance X/H is taken from the equation 



* E-mails: v.escalante@crya.unam.mx; 
chris.morisset@gmail.com; georgiev@astro.unam.mx 



/(A) /e(A) dV (e(A)) X 
/(H/3) Je(Hp)dV (e(H/3)) H 



(1) 
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where I (X) / 1 (H/3) is the ratio of intensities corrected for 
reddening, e(A) and e(H/3) are the emissivities of a line emit- 
ted by X and of H/3 respectively, both are integrated over 
some volume fraction of the nebula covered by the aperture, 
and (e) means some kind of average of the emissivities over 
that volume fraction. Averaged emissivities over the volume 
within the solid angle of the aperture imply the adoption 
of average values for the temperature, density and filling 
factors. 

An alternative procedure to determine chemical abun- 
dances in a nebula is to propose mutually consistent mod- 
els of the nebula and the atmosphere of the exciting star 
or stars. The stellar parameters, chemical abundances and 
density profile of the nebula are inputs to a photoionization 
model, and their values are fixed by matching the predicted 
line intensities to the observations. The variation of stellar 
and nebular parameters of the model can be constrained 
by direct observations of the star and nebula to reduce the 
degeneracy of different param eters giving similar matche s 
to the observed line intensities. iMorisset fc Georgiev! 4.2009) 
(hereafter MG) showed that a consistent stellar energy dis- 
tribution (SED) of the central star of IC 418 and a photoion- 
ization model give an impressive agreement of predictions 
and measured intensities of CELs and many recombina- 
tion lines within observational errors. However most dipole- 
allowed lines in the SWBH survey were not included in the 
model by MG, and some of the ones that were included show 
observed intensities much larger than their predicted values 
because the MG calculations did not include fluorescence 
excitation. 

In this work we add fluorescence excitation of permit- 
ted lines to a photoionization model and a SED similar to 
the ones used by MG, and compare the results with all lines 
observed by SWBH that have probable identifications as 
dipole-allowed lines of Cn, Nl, Nil, O I, and On and op- 
tical forbidden lines of N I and O I. Continuum fluorescence 
of starlight was suggested as a mechanism to e xcite permit- 
ted lines in i onized nebulae by ISeatonl dl968h and iGrandil 
l|l975l . ll976l '). lEscalante and Morissetl (|2005l ) (hereafter EM) 
showed that fluorescence in Orion dominates the intensities 
of N n permitted lines. The low excitation level of IC 418 and 
the resemblance of its ionization degree to that of the Orion 
nebul a noted by iTorres-Peimbert. Peimbert fc Daltabuitl 
( 1980) suggest that a similar fluorescence mechanism may 
occur in both objects. 



2 LINE EMISSIVITIES 

In our previous work (EM) we found that the rate of contin- 
uum fluorescence excitation typically depends on the absorp- 
tion of photons in several resonant transitions, also called 
pumping transitions, to quantum states of moderate excita- 
tion. Excitation by recombination on the other hand involves 
the capture of free electrons into a much larger number of 
excited states. In the fluorescence and recombination mech- 
anisms subsequent transitions to lower states produce the 
optical lines observed in ionized regions although many lines 
in the UV and IR are also emitted. 

Our calculation of quantum state populations is based 
on the cascade matrix formalism described in detail by EM. 
In that formalism excited states are populated by recombi- 



nations and by transitions from the ground and metastable 
states. The contribution of recombinations to the population 
of a state j is given by the effective recombination coeffi- 
cient of the state, a^ ff . The effective fluorescence coefficient 
defined by EM, (3gf, gives the contribution of pumping tran- 
sitions from the ground and metastable states to excited 
states that produce decays ending in state j. Both coeffi- 
cients can be calculated in terms of the cascade matrix of 
ISeatonl l|l959f ) . The emissivity of a line for a transition from 
state j to i with a branching ratio Pji and frequency v at a 
point in the nebula is 

e 3l = huP^n+af + n 3 Pfj) (2) 

9 

where n e is the electron density, n + is the density of the 
ion before recombination, and n g is the population density 
of the ground or metastable state of the recombined ion or 
atom. The branching ratio Pji = Aji/Aj, where Aji is the 
Einstein coefficient for transition from states j to i, and Aj 
is the total decay rate of the upper state j. The line effective 
recombination coefficient is defined as: 

a ji = p n a T ( 3 ) 
It is assumed that the ion density n + can be calculated 
from the ionization balance condition, and the ground and 
metastable state populations can be calculated from the bal- 
ance of collisional and radiative transitions with the aid of 
a nebular photoionization model. 

The cascade matrix formalism assumes that every tran- 
sition from the continuum or from the ground or metastable 
state to a higher-energy state is followed by transitions to 
lower-energy states, and that metastable states are not sig- 
nificantly affected by those transitions. Those are valid as- 
sumptions at the temperatures and densities of a photoion- 
ized nebula because collisional excitation and deexcitation 
with electrons control the population of metastable states, 
but fluorescence can compete with collisions with electrons 
and atoms in the cold er interfa c e of t he neutral and ion- 
ized part of a nebula. iBautistal ([l999) proposed that the 
forbidden lines of atomic N are excited by fluorescence and 
collisions with electrons and atoms in Orion. Because the 
collisional excitation of metastable states violates the as- 
sumption of downward transitions in the calculation of the 
cascade matrix, the population densities of the metastable 
and ground states, rij and rife, must be solved for from equa- 
tions of the type: 

ns + = J2 Uk ( Ak > + Ck i + ft*) (4) 

where Aji — if state i is energetically higher than j, and 
Cji is the collisional excitation or deexcitation rate from 
states j to i. The summations comprise all the metastable 
and ground states, except state j. The system of equa- 
tions @ is non-linear because the effective fluorescence 
rates depend on the optical depth of transitions from the 
ground and metastable states to excited states (see equa- 
tions [2] through [7] of EM). We found that the sys- 
tem (|4]) can be solved iteratively by updating the optical 
depths across the nebula, and then using matrix inversion 
to solve for the populations rij at each point in the neb- 
ula in each iterat i on. T he five-state inversion formulae by 
iKafatos fc Lynch! |l99d) can be easily modified to include 
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the fluorescence rates in equation Q. The emissivity of a 
forbidden line is then 

e C ji = hvAjiUj (5) 

Effective recombination coefficients depend on tempera- 
ture and more weakly on density. They have been tabulated 
with their radiative and dielectronic parts for spectroscopic 
terms of astrophysically important ions in the literature. We 
assume that af 1 for an individual fine structure state of a 
given term is proportional to the relative weight of the state 
(2J+1) / JJ.(2Jj+l). This is a good approximation for direct 
recombinations from the continuum. Nevertheless it is worth 
noticing that o^ ff adds the contribution of a large number 
of cascading bound-bound transitions after the free-bound 
transition, and probabilities of fine structure transitions are 
not proportional to relative weights. To date only the calcu- 
lation of effective recombination coefficients of N + by FSL 
has considered the fine structure splitting of quantum states. 

In the calculation of /3gf we split quantum states in their 
fine structure components either assuming the LS or LK an- 
gular momentum coupling schemes as appropriate or using 
published intermediate coupling calculations when available. 
The fine structure of states can be important in fluorescence 
calculations because quantum selection rules can selectively 
pump certain states. 

The calculation of /3gj involves a large number of quan- 
tum states that can be reached by downward transitions 
following an absorption from the ground or a metastable 
state. Metastable states and some core-excited states of 
the form 2s2p m+1 can build up a significant population 
under low density conditions because of their low tran- 
sition probabilities to the ground state. The population 
densities, n g , in equation (0 include all the fine struc- 
ture states of the ground and metastable terms. Optical 
lines involving core-excited configurations have often been 
attributed to low temperatu re dielectronic recombination 
|Nussbaumer fc Storev|[l98ll '). but radiative recombination 
or fluorescence followed by two-electron transitions can also 
produce excited-core configurations. EM found that it is im- 
portant to take into account spin forbidden transitions and 
core excited configurations in the calculation of /3gJ. 



3 ATOMIC DATA 

We used the compilation in the NIST Atom ic Spectra 
Database (Ralchcn ko. Kramida fc Read er 2008) to classify 
quantum states and store the observed energies and proba- 
bilities for some transitions. However we relied substantially 
on more recent and complete calculations of probabilities 
from other authors. Systematic variations among different 
atomic databases have small effects in the cascade matrix 
in most cases because the calculation depends on branching 
ratios. 

3.1 C+ 

This ion has the simplest structure. The ground core term 
produces a doublet manifold of excited states 2s 2 ( 2 S) nl 2 L. 
Consequently free-bound transitions tend to concentrate in 
fewer states, and give rise to more intense lines from states 
with higher quantum numbers. The excited cores 2s2p, and 



2p 2 produce quartet and doublet manifolds. We used the 
populations of the fine structure states of terms 2s 2 2p 2 P° 
and 2s2p 2 4 P as the variables n g in equation @ for this 
ion. Some excited core configurations, like 2s2p( 3 P°) 3s and 
2s2p( 3 P°) 3p, were also detected in IC 4f8 by SWBH. 

We used te r m-ave raged transition probabilities calcu- 
lated by iNaharl (|l995h with the R-matrix approach sup- 
plemented with spin forbidden transitions from the NIST 
databas e. The latter compiles t he Op acity Project calcula- 
tions bv lYan. Taylor fe Seatonl dl987T). and some int e rmed i- 
ate coupling calculations by iNussbaumer fe Storevl <|l98lf ). 
We split transition probabilities amo ng fine struc ture com- 
ponents assuming pure LS coupling l|Allenlll973l ) since no 
significant departures from th at coupling have been reported 
l|Wiese. Fuhr fc DetersHl988h . We used the effective recom- 
binati on coefficients calculated bv lDaval. Storey fc Kisieliud 
|2000l ). 

3.2 N+ 

The ground-core term of this ion produces singlet and triplet 
manifolds of excited states. The fine-structure populations 
of the terms 2s 2 2p 2 3 P, X D and 1 S, were used as the pop- 
ulations n g in equation (J2J. Terms with the excited cores 
2s2p 2 4 P and 2 D were detected in IC 418 by SWBH. Those 
cores produce singlet, triplet and quintet manifolds. The 
quintets are not reachable by transitions from the continuum 
if LS coupling selection rules prevail. Therefore the lines ob- 
served by SWBH between AA 5526.2 and 5551.9, which are 
probably identified with the Nil multiplet 2s2p 2 ( 4 P)3s 5 P- 
2s2p 2 ( 4 P)3p 5 D°, may be produced by dielectronic recombi- 
nation out of LS coupling followed by stabilizing transitions 
to the 3p D° term. Recombination coefficients for those 
transitions have not been calculated yet to our knowledge. 

The atomic data to calculate this spectrum was de- 
scribed in detail by EM. Different calculations of transition 
probabilities — with due account of the fine structure split- 
ting and breakdown of t he LS coupling scheme in 4f an d 5g 
states — generally agree l|Lavm. Olalla fc Martmll2000l ). 

Most calculations of effective recombination coefficients 
give rate s for terms in some coupl ing scheme, usually LS 
coupling. lEscalante fc Victor! (|l990T i (hereafter EV) calcu- 
lated ctji for the f states assuming an LK coupling scheme. 
In that scheme, the core total spin, S p = 1/2 is added to the 
total angular momentum V to produce an intermediate an- 
gular momentum K = \S P — L'\ . . . \S P + L'\ with values 3/2, 
5/2, 7/2 and 9/2 for the f states. The valence electron spin 
is added to K to give a doublet of total angular momentum 
J = K ± 1/2 for each LK term. IVictor fc Escalante I (|l988l ') 
and EV give transition probabilities and a^f averaged over 
the K states of each mf L' term. To get the transition prob- 
ability to a given K term, the /^-averaged probability of 
the transition rod ' 3 L — > mf L' must be multiplied by the 
relative weight (2K + 1)/(2S P + 1)(2Z/ + 1) of the term. The 
probability of the inverse transition mf L' — s> nd 1,3 L is the 
same for all the K states of the mf L' term. We split the 
3d 1,3 L-4f L' mul tiplets into fine structure compo nents with 
formulae given bv lEscalante fc G6ngora-T.I l|l990f ). Effective 
recombination coefficients of LK terms were multiplied by 
the relative weight, (2 J+1)/(2S I P + 1)(2L' + 1), to get the fine 
structure ocf 1 ■ All of the above calculations to obtain the re- 
combination and transition rates to individual fine structure 
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Table 1. Comparison of effective recombination coefficients and A values for selected Nil lines at 10 4 K , (f) Branching r atio o f the 
line or multiplet from EV , (2) statistical weight of term or state, (3) FSL, (4) EV, (5) iMar et alj ||2000| 1, (6) IShen et al.l <201C|) , (7) 
IVictor fc Escalante I il988h . 



Transition 


A (A) 


p.. 
(1) 


q 
(2) 


a cff 
(3) 

n e = 10 3 cm 3 


(10 -15 cm 3 s~ 
(3) 

10 4 cm" 3 


■1) 

(4) 


(5) 


Aji(10 8 s" 
(6) 


- 1 ) 

(7) 


3p 


3 D-3d 3 F° 


5006 


0.942t 


21 






238.00+ 








3p 


3 D 3 -3d 3 F° 


5005.15 


1.000 


9 


77.20 


105.00 


108.28 


— 


— 


1.246 


3d 


3 D3-3d 3 F§ 


oOzo.oo 


0.092 


7 


8.69 


7.92 


7.78 


0.118 




0.107 


3p 


3 D 2 -3d 3 F° 


5001.47 


0.907 


7 


85.60 


78.00 


76.36 







1.050 


3p 


3 D 2 -3d 3 F° 


5016.30 


0.140 


5 


9.60 


7.43 


8.42 


0.186 


— 


0.162 


3p 


3 Di-3d 3 F° 


5001.13 


0.844 


5 


58.20 


45.00 


50.75 


— 


— 


0.976 


3d 


3 P°-4f D 


4433 


0.416 1 " 


20 






27.20t 








3d 


3 Pg-4f D[3/2]i§ 


4433.48* 


0.461 


3 








4.52 


— 


1.300 


1.070 




3po ac T-^rc /n] 
Jr^ — 41 U[D/ Z\2 


4442.02 


0.258 


5 


5.17 


5.35 


4.22 


0.695 


0.824 


0.578 


3d 


3 P§-4f D[5/2] 3 


4432.74* 


0.830 


7 


18.40 


19.00 


18.98 


— 


1.790 


1.926 


3d 


1 P°-4f D 


4690 


0.129t 


20 






8.46t 








3d 


1 P°-4f D[3/2] 2 


4678.14* 


0.317 


5 






5.19 




1.030 


0.717 


3d 


ipj-ytf D[5/2] 2 


4694.64 


0.481 


5 


6.18 


6.40 


7.88 


0.607 


0.753 


1.076 


•~±A 

OQ 


lno /\f -p 
u — 41 r 


4174 


0.182t 


28 






18.30t 








3d 


!D°-4f P[7/2] 3 


4171.60 


0.508 


7 


16.50 


12.10 


12.76 


0.448 


0.544 


1.181 


3d 


!D°-4f F[5/2] 3 


4176.16* 


0.380 


7 


13.20 


9.18 


9.54 


1.130 


1.280 


0.886 


•~±A 
OQ 


3no /if "p 
u — 41 r 


4242 


0.554t 


28 






55.60 1 " 








3d 


3 D°-4f F[5/2] a 


4236.93* 


0.747 


5 


21.00 


15.00 


13.39 




1.850 


1.758 


3d 


3 D§-4f F[5/2] 2 


4246.86? 


0.004 


5 








0.07 





0.007 


0.009 


3d 


3 D°-4f F[5/2] 3 


4241.76* 


0.456 


7 


20.60 


14.40 


11.43 


— 


0.728 


0.886 


3d 


3 D§-4f F[7/2] 4 


4241.79* 


0.889 


9 


36.10 


26.80 


28.69 




1.960 


2.090 


3d 


3 D°-4f F[7/2] 3 


4237.05? 


0.343 


9 


38.58 


6.29 


8.60 




0.986 


0.797 


■'J. A 
OQ 


3t?0 AC T? 

r —41 r 


4089 


0.069t 


28 






6.92t 








3d 


3 F°-4f F[7/2] 4 


4095.90 


0.104 


9 


2.77 


2.05 


3.35 





0.144 


0.244 


3d 


3 F°-4f F[7/2] 3 


4096.57 


0.004 


7 


0.103 


0.075 


0.10 








0.009 


3d 


3 F§-4f F[7/2] 4 


4082.27 


0.007 


9 


6.86 


5.08 


0.22 


0.335 


0.485 


0.016 


3d 


3 F°-4f F[7/2] 3 


4082.89 


0.043 


7 


0.41 


0.30 


1.08 




0.011 


0.100 


3d 


3 Fg-4f F[7/2] 3 


4073.05 


0.004 


7 


7.90 


5.81 


0.096 


0.499 


0.631 


0.0089 


3d 


1 F°-4f G 


4543 


0.234t 


32 






28.30t 








3d 


*F°-4f G[9/2] 4 


4530.41* 


0.535 


9 


15.30 


19.90 


18.18 


1.450 


1.490 


1.242 


3d 


*F°-4f G[7/2] 4 


4552.52* 


0.424 


9 


9.10 


8.53 


14.41 


0.611 


0.650 


0.993 


3d 


3 F°-4f G 


4039 


0.765 1 " 


32 






92.40+ 








3d 


3 F°-4f G[9/2] 4 


4026.08 


0.436 


9 


7.15 


9.34 


14.81 


0.672 


0.794 


1.013 


3d 


3 F°-4f G[9/2] 5 


4041.31* 


0.999 


11 


28.60 


37.90 


41.48 


2.080 


2.780 


2.430 


3d 


3 F°-4f G[7/2] 4 


4043.53* 


0.540 


9 


18.00 


16.90 


18.34 


1.250 


1.130 


1.266 


3d 


3 F°-4f G[7/2] 3 


4035.08* 


0.917 


7 


19.10 


17.90 


24.24 


1.300 


1.650 


2.232 


3d 


3 F°-4f G[7/2] 3 


4058.16? 


0.0012 


7 


0.080 


0.074 


0.033 






0.0031 



§ The notation in LK coupling is 4f 
* Observed in IC 418. 
? Uncertain identification in IC 418. 
t a ' 1 averaged over K values. 



states assume that the radial part of the dipole transition 
element of an LS or LK term does not change significantly 
among different J states of a term. 

Recently FSL calculated effective recombination coeffi- 
cients of Nil fine structure states in intermediate coupling 
and including the effects of dielectronic recombination. Their 



rates of certain lines show a hitherto unknown strong depen- 
dence on electron density due to the variation of the popu- 
lations of the parent ion ground states with density. Table [1] 
shows a comparison of effective recombination coefficients at 
10 4 K and two electron densities for some of the most intense 
lines excited by recombination at that temperature. 
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There is general agreement between the new effective 
recombination coefficients of FSL and those of EV at the 
electron densities of n e — 10 3 cm -3 and 10 4 cm -3 , but there 
are significant differences in some transitions, most notably 
in multiplet 3d 3 F-4f F[7/2]. Since there is agreement be- 
tween the two calculations in some transitions of that term, 
the differences in the other transitions with the same upper 
state appear to arise from the branching ratio of the transi- 
tion, Pji , or more specifically from the transition probability. 
For the transiti ons with the largest disagre ements in a|f , the 
probabilities of lVictor fc Escalante I (|1988T) are much s maller 
than the experimental measure ments by iMar et al.l (|2000| ) 
and relativistic calculations by IShen. Yuan fc Liul |2010l ). 
Those discrepancies may point to important relativistic ef- 
fects in some transitions th at were not taken into accou nt in 
the term-averaged values of IVictor fc Escalante I (1 19881 ) and 
the way in which we split the probabilities among the fine 
structure components. Unfortunately those transitions were 
not detected in IC 418 by SWBH, although they should be 
just below their detection limit if the larger branching ratios 
are correct. 



3.3 0+ 

This ion has a ground configuration 2s 2 2p 3 . The populations 
of the fine-structure states of its terms 4 S°, 2 D° and 2 P° 
were used in equation ([2]) as the n g variables. Those terms 
are connected to quartet and doublet manifolds of excited 
states with several cores. SWBH possibly identified 85 tran- 
sitions with the 2s 2 2p 2 3 P ground core, 32 transitions with 
the 2s 2 2p 2 X D excited core, 12 transitions involving both 
cores, and a few lines with the 2s2p 3 5 S° exci ted cor e. 

We used the transition probabilities of iNaharl (I2009T I 
supplemented with th e intermed i ate co upling calculations 
for the 3d-4f array by iLiu et al.1 (|l995h . who found signif- 
icant effects of LS coupling breakdown in 4f and some 3d 
terms. Terms with I > 3 are better described in LK cou- 
pling (Wenakerl ll990l). Effective recombination coefficients 
were take n from ILiu et al.l |l995) for the 3d and 4f terms, 
and from IStorevI (|l994i r for other states. The effective re- 



states were taken from Berrington fc Burke! (|l98ll ) and 
iBell. Berringtpn fc Thomas (119981) for electron colli- 



and 



Krems. Jamieson fc Dalgarnol (|2006l ) and 



lAbrahamsson. Krems fc Dalgarnol ( 2007 ) for collisions with 



combination coefficients of lNussbaumer fc Storevl (|l984T ) for 
terms with the 2s 2 2p 2 X D core involve only the dielectronic 
contribution. Transitions between states with the 2s 2 2p 2 3 P 
core and states with the 2s 2 2p 2 X D core, some observed with 
non-negligible intensity in IC 418, may add a contribution 
of radiative recombination to the terms with the 2s 2 2p 2 *D 
core at lower temperatures. 



3.4 N and O 

For nitrogen we used the NIST compilation of transition 
probabilities mentioned at the beginning of this section. 
For oxygen we used the model potential calculations of 
lEscalante fc Victoij (|l994T i complemented with the NIST 
compilation for spin forbidden and two-electron transitions. 
Since these atoms are abundant in the neutral zone, 
the main transitions pumping the lines by fluorescence 
will be below the Lyman limit. E ffective recombina- 
tion coeffici ents were taken from lEscalante fc Victor! 
1 19921) and iPeauignot. Petitiean fc Boissonl (Il99ll ). Ef- 
fective collision strengths for excitation of O metastable 



H atoms. Effective collision strengths for excitati on of 
N me tastable states by electrons were taken from iTavall 
(2006). Transition pro babilities for the N and O metast able 
states were taken from lFroese Fischer fc Tachievl |2004h . 

Atomic N has the same atomic configurations as + and 
the populations of the five states of the metastable terms 
4 S°, 2 D° and 2 P° were found from equation ((4}. Atomic O 
has a ground configuration 2s 2 2p 4 with terms 3 P, *D and 1 S, 
and its excited states in the ground core configuration form 
quintet and triplet manifolds. There is a singlet manifold 
with a 2s 2 2p 3 ( 2 D°) excited core. SWBH detected 10 possible 
lines with that core, but there are no calculations of effective 
recombination coefficients for them. Recombination from the 
ground states of + or fluorescence from the 3 P states can 
affect the populations of the metastable singlet states D 
and S through spin-forbidden transitions from the excited 
triplet states. Those transitions will be important in the cas- 
cade only if the stronger resonant trans itions of the triplets 
are o ptically thick, the so called case B (jEscalante fc Victor] 
1992), and the temperature is low enough to prevent colli- 
sions with electrons and H atoms from dominating the pop- 
ulations of the O metastable states. 



4 MODEL CALCULATIONS 

The efficiency of fluorescence to enhance dipole-allowed lines 
depends on the transfer of starlight, the intensity of the dif- 
fuse field and the ion distribution throughout the nebula. 
Therefore a detailed SED and model of the nebula are re- 
quired to predict adequately the intensity of permitted lines. 



4.1 Stellar energy distribution 

To evaluate the amount of emitted energy better and in- 
crease the resolution of the SED near the wavelengths of 
the most important pumping transitions, we calculated a 
new model of the central star atmosphere using the CMF- 
GEN code (Hillier and Miller, 1998) with a more complex 
atomic structure than the one used in MG as shown in ta- 
ble [5] The increased number of lines in the far UV region 
of the spectrum changed the ionization balance of the neb- 
ula. Following the same procedure used by MG, we adjusted 
the stellar atmosphere model, increasing its temperature to 
39000 K and keeping the other parameters unchanged. The 
new self-consistent stellar and nebular models reproduce the 
intensities of the emission lines common to both models with 
a quality factor k(O) (see definition in equation (6) of MG) 
of the same order of magnitude as the one obtained with the 
previous SED. The new high-dispersion (R ~ 100000) SED 
was corrected for a stellar rotation velocity of 30kms _1 as 
derived in MG. 

The Doppler shift of the resonant transition frequencies 
due to internal motions of the nebular gas must be taken 
into consideration with a SED of that resolution as discussed 
in section 15.61 We adopted a velocity field V oc R 4 with a 
maximum velocity V m = 40kms~ 1 at the outer edge of the 
nebula. MG adopted a similar velocity law to reproduce the 
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Figure 1. Continuum-normalized SED with identifications of stellar lines. Unlabelled lower vertical lines (blue) show the wavelengths 
of the transitions in Nil multiplet 2p 2 3 P-3d 3 D°AA 533.5f-533.88, and Nil line 2p 2 1 D 2 -3d !D°A 582.156 in the rest frame of the 
photosphere. 



profiles of the H/3, [Nil] and [0 ill] lines (see appendix B of 
MG). The effects of other velocity fields and of turbulence 
on the calculated line intensities are discussed in section [5U1 

We did not find a systematic decrease of fluorescence 
excitation of the CNO lines in the SWBH survey with a 
higher-dispersion SED as reported for the Balmer lines by 
iLuridiana et al.l |2009f ) . The probable reason is the fact that 
excited states in multielectron species have several resonant 
transitions distributed over a range of wavelengths. The flu- 
orescence excitation to those states can increase or decrease 
with a high-dispersion SED depending on the coincidence 
of those resonant transitions with the emission peaks or ab- 
sorption troughs of the stellar wind profile in the rest frame 
of the absorbing ion or atom. The resolution of the adopted 
SED is larger than the width of the pumping lines in the 
gas. Therefore we do not expect any remaining bias in the 
predicted intensities with the adopted SED. 

Fluorescence excitation depends on a large number of 
resonant transitions, but there are a few transitions espe- 
cially important in the populations of certain states. Fig. [T] 
shows two sections of the SED with examples of important 
transitions in the excitation of N II lines by fluorescence. The 
low abundance of N + in the stellar atmosphere produces lit- 
tle depression of the continuum at the Nil pumping tran- 
sitions in the gas, but other species do produce important 
absorption features in the SED as shown in Fig.[Tl 

4.2 Photoionization model 

The photoionization model calculated by MG took into ac- 
count the small ellipticity of IC 418 by using the pseudo-3D 



code CLOUDY_3D (|MorissetJl200fJ ). which allows the calcula- 
tion of nebular models with complex geometries. The near 
spherical symmetry of IC 418 and the small angle of the 
aperture used by SBW allowed us to approximate the neb- 
ula with a ID spherical model to speed up the calculation of 
thousands of transitions in this work. We extracted contin- 
uum opacities, ground and metastable state populations, ion 
and electron densities, and electron temper atures from ver- 
sion c08.00 of CLOUDY l|Ferland et all ll99Sf ) with the SED 
described above and a density distribution similar to the 
one used by MG. Emissivities calculated with equations ([2]) 
and ((5]) are then integrated along a beam across the nebula 
to simulate the size and position of the aperture used by 
SBW. Images of the slit position on the nebula are given in 
MG. 

The stellar light is attenuated by geometrical dilu- 
tion and continuum opacity as given by CLOUDY. The self- 
absorption within the resonant transition is taken into ac- 
count with the 'pumping probability' denned by iFerlandl 
(1992). The calculation of ffi;,- uses th e escape probability 
formalism of iHummer fc Storevl |l992) to approximate the 
propagation of photons in the resonant transitions with an 
overlapping continuum in the nebula. This in effect takes 
into account the variation of optical depth in a medium par- 
tially thick in the resonant t ransit ions, a situation named 
'case D' by ILuridiana et all (2009). Additional details can 
be found in EM. 

Although MG used some nebular lines to adjust the 
stellar atmosphere and nebular model parameters, most of 
their calculated and all of the calculated line intensities in 
the present work are ab initio predictions of the SED and 
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Table 2. Ions used in the stellar model 



Table 3. Element abundances in the nebular model 



Ions 


Levels 


Superlevels 


Hi 


30 


20 


He i 


69 


45 


Hen 


30 


30 


Cn 


92 


40 


Cm 


243 


99 


Civ 


64 


59 


Nn 


85 


45 


Nm 


70 


34 


Niv 


76 


14 


Nv 


67 


45 


On 


123 


54 


Oin 


170 


88 


Oiv 


78 


38 


Ov 


152 


75 


Ovi 


13 


13 


Sim 


34 


20 


Si iv 


33 


22 


Piv 


178 


36 


Pv 


62 


16 


Sin 


28 


13 


Siv 


112 


51 


Sv 


98 


31 


Svi 


58 


28 


Nc ii 


18 


14 


Nem 


71 


23 


Neiv 


52 


17 


Arm 


36 


10 


Ar iv 


105 


31 


Ar v 


99 


38 


Fein 


1433 


104 


Feiv 


1000 


100 


Fev 


1000 


139 


Fe VI 


433 


44 


Fe VII 


153 


29 



nebular model. The abundances adjusted by MG with this 
procedure and used in this work are given in table The 
density profile in this work is 



n(R < R in ) 
n(R 3 > R > R in ) 







(6) 



n + n\ exp \(R — i?i) 2 /o"i] + 

n 2 exp [(R - R2) 2 Vol] 

where R is the distance to the central star in cm. Table 0] 
gives the adopted values for parameters no. ..2, Rin, Ri, R2, 
<7i and 02. Parameters R\, R2 have intermediate values be- 
tween the long and short axes in the ellipsoidal model of 
MG. This density profile was chosen to match the HST im- 
ages of the intense nebular lines as discussed in detail by 
MG (see their figure 5). With few exceptions, the predicted 
line intensities in our ID model vary by less than 30 per 
cent when no. ..2 are varied simultaneously by 20 per cent or 
when Ri and R2 are varied simultaneously by 20 per cent. 

4.3 The photodissociation region 

It has long been known that IC 418 has a PhotoDissoci- 
ation Region (PDR) extending beyond the opti cal nebula 
(|Cohen fe Barlow! Il974l ; I Taylor fe Pottaschlll987T) with low 
ioniz a tion and atomic compo n ents ( Monk. Barlow fe Clegd 



Il99d ; iTavlor. Gussie fe Gossl Il989l : 
© 0000 RAS, MNRAS OOO.rrH20l 



Meixner etafl 1996) 



Element 


logX/H 


He 


—0.92 


c 


-3.10 


N 


-4.00 


O 


-3.40 


Ne 


-4.00 


Mg 


-4.95 


Si 


-4.90 


S 


-5.35 


CI 


-6.90 


Ar 


-5.80 


Fe 


-7.40 



Table 4. Nebular density profile parameters 



Parameter Unit 



Value 



no 
m 
■n-2 
"3 

log R in 
log Ri 
logi?2 
log #3 
log m 
log cr 2 



cm 
cm 
cm 
cm 
cm 
cm 
cm 
cm 
cm 
cm 



2850 
5260 
9500 
20500 
16.06 
16.28 
17.02 
17.10 
15.96 
16.46 



ITavlor et all (|l989h measured an H mass of 0.35 ± 0.05 M© 
over a radius > 90 arcsec from the center of the nebula. 
Thus far there have been no detections of molecular emis- 
sion in IC 418 llHuggins et al.lll996l: iDaval fe Biegindll996l : 
iHora et all Il999l ). (Liu et al.l (|2001aJ ) obtained an average 
density for the PDR of 10 5 ' 5 cm~ 3 from observations of the 
[Cn] 158- /im and [Ol] 63- um and 146- um lines. That den- 
sity is much higher than the one needed to explain the radio 
observations and suggests the existence of a narrow dense 
shell just outside the ionization front or clumps in the PDR. 

The CLOUDY nebular model can be readily extended to 
the PDR by adding an outer neutral atomic envelope to the 
ionized region starting at a distance of the central star R3. 
We chose the following density profile for the PDR: 



n PDR (fl > R 3 ) = n 3 (T e /10 4 K)" 



(7) 



where T e is the electron temperature in kelvins, and n 3 and 
R3 are given in table [4] This density profile is constrained 
by the observed fluxes of the FIR [C n] and [0 1] lines. The 
temperature dependence in Eq. © allows for a smooth tran- 
sition between the ionized and neutral material by increasing 
the density as the temperature drops. The extent of the neu- 
tral region needed to match the observed intensities of the 
FIR 1 and C 11 lines is just a few percent of the radius of 
the ionized region. A more realistic model of the ionized and 
molecular interface needs a hydrodynamic model, which is 
beyond the scope of this paper. 
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5 RESULTS 

5.1 Selection of observations 

Proper line identification of weak spectral lines is a funda- 
mental problem in the analysis of deep nebular spectroscopy. 
SWBH presented the line identification code emili to auto- 
mate and standarise an otherwise tedious labour of iden- 
tifying several hundreds of spectral lines. SBW selected a 
subset of emission line intensities from the survey of SWBH 
based on the likelihood of a correct, unique identification 
as given by emili and the availability of published effective 
recombination coefficients to determine abundances. Some 
lines in SWBH were not taken into account by SBW be- 
cause those lines were probably blended with other emission 
lines or emili gave more than one plausible identification for 
them. 

We have reexamined the SWBH line list and added 
many measured lines of Cn, Nl, Nil, Ol and On to the 
selection of SBW in order to compare them with our cal- 
culations either because we could estimate their effective 
recombination coefficients reasonably well or because they 
were the most probable contributor to a blend. Some cases 
for individual species are discussed in the following sections. 

Published calculations often tabulate the line effec- 
tive recombination coefficient in equation © for the most 
intense lines only. Coefficients for weaker lines with the 
same upper state j can be readily estimated by dividing 
by the branching ratio Pji of the tabulated transition and 
multiplying by the corresponding Pji of the desired line 
as in equation In order to do this consistently, it is 
necessary that authors publish the transition probability 
data or at least the branching ratios used in their calcula- 
tions llNussbaumer fc Storevl 19841; lEscalante fc Victor! 1990l; 
Pequignot et all Il99ll ; lEscalante fc Victorl Il992l ; iLiu et ail 
19951 ). 

The issue of multiple identifications is more compli- 
cated. A line with multiple identifications may be either a 
real blend or some of the proposed identifications may not be 
real. Our calculations show that in many cases a transition in 
a possible blend contributes insignificantly to the measured 
intensity. When a transition in a blend of two transitions 
of the same species contributes more than 10 per cent of 
the total intensity, we have added the calculated intensities 
of the transitions to compare them with the measurements 
by SWBH. Blends of emissions of the same species are thus 
compared with our computed spectra without the need to 
assume theoretical intensities to 'deblend' the observed lines. 
Blends of emissions of different species are compared directly 
with the predicted intensities. 



5.2 Comparison with observed intensities 

In Fig. [2]to[S]we compare predicted intensities with observed 
intensities in a scale of /(H/3) = 10 4 . Tables [51 to 1 1 1 1 give the 
predicted intensities of dipole-allowed lines with intensities 
higher than 8 x 10~ 6 of the H/3 strength, which is close to 
the lower limit of detection in SWBH. Tables [lOl and [TT1 also 
include predicted intensities of the optical forbidden lines of 
N I and O I. Some of the wavelengths given in the second col- 
umn of the table are calculated from energy levels and have 
limited accuracy. The fraction of calculated intensity due to 



recombination is given in the third column of the tables. 
The calculated and observed intensities in the fourth and 
fifth columns are normalized to /(H/3) = 10 4 . The bottom 
line in each table gives the average of the ratio of the inten- 
sity due to recombination to the total calculated intensity 
and the average of the ratio of the calculated intensity to 
the observed intensity for each species. We have found some 
undetected lines with calculated intensities above the detec- 
tion limit in SWBH. Some of them are blended with intense 
nebular lines, and a few were measured, but apparently not 
identified by emili. The observed intensity of lines not se- 
lected by SBW are marked with '?' in the sixth column of 
the tables. 

The uncertainty of each line intensity depends on sev- 
eral factors like the signal-to-noise ratio (S/N), its relative 
width with respect to the instrument resolution and whether 
the line may be blended with other lines. In order to make 
meaningful comparisons with observations, MG estimated a 
variable uncertainty of the line fluxes in SBW between 10 
per cent for the brighter lines to 30 per cent for the weaker 
lines. SBW suggest an uncertainty of 20 per cent for lines 
with a S/N > 20. Most of the lines analysed in this work 
are weaker than 0.01 times the intensity of H/3, and have 
a S/N > 20. Therefore we will adopt a minimum general 
uncertainty of 20 per cent for all dipole-allowed lines in the 
SWBH survey. 

Our calculations tend to underestimate mo st of the lines 
with d ubiou s identification according t o emili . iRola fc Pelal 
i| 19941) and IWesson. Stock fc Sciclunal l|2012t ) have shown 
that measured intensities of lines with S /N < 6 are strongly 
biased upwardly. We only have two lines with an underesti- 
mated predicted intensity and a measured S/N < 7, which 
corresponds to a true S/N between 6 and 8 at the 1-cr con- 
fidenc e level according to the calculations of IRola fc Pelal 
1 19941) . In most of our cases it is more likely that the un- 
derestimated predictions indicate incorrect identifications or 
blends with other lines rather than uncertainties in the mea- 
sured values. We have not included in the tables and figures 
features with a dubious identification and a predicted inten- 
sity less than 0.1 of the measured value. 

The excitation mechanisms of a spectrum can be diag- 
nosed by contrasting the intensities of lines yie l ded by differ - 
ent atomic configurations as done by iGrandil (|1975l . |l976j) . 
Lines from s, p and most d states with the same spin mul- 
tiplicity as that of the ground state are enhanced in vary- 
ing degrees by fluorescence relative to recombination lines if 
the recombined ion concentration is sufficiently high. Tran- 
sitions from d states that are not connected to the ground 
or a metastable state by a resonant transition and transi- 
tions from f and g states are excited mostly or totally by 
recombination. 

SBW found that some dipole-allowed lines of a species 
have profiles similar to those of the forbidden lines of the 
next higher ionization stage of the species, (e.g. some O I line 
profiles are similar to [O n] profiles), indicating that they are 
produced in the same part of the nebula by recombination 
of the more ionized stage. Likewise they found that other 
dipole-allowed lines of the same species have profiles similar 
to those of the forbidden lines of the same ionization stage, 
indicating that they are produced by fluorescence excitation 
of that species. 

SBW also found a correlation between the line width of 



© 0000 RAS, MNRAS OOO.ni[20l 



Fluorescence and recombination in ICJf.18 9 



the species producing the line and the ionization potential 
of the species, which indicates a relationship between the 
line width and the ionization stratification in the nebula. 
Recombination lines of species with lower ionization poten- 
tial generally have broader lines and their widths are similar 
to those of the forbidden lines of the next higher ionization 
stage, indicating again a common place of origin within the 
nebula. 

SBW proposed that line profile and line width similar- 
ities can be used to discriminate the excitation mechanism 
of the lines. The contributions of fluorescence and recombi- 
nation, however, are comparable in many lines and produce 
many exceptions to the discriminating criteria proposed by 
SBW as discussed for particular cases below. Therefore a 
more quantitative analysis is needed to determine the exci- 
tation mechanism of a line. 

5.3 Cn 

The agreement of predicted recombination intensities within 
20 per cent of observed values in the SBW selection of 9 lines 
or blends of lines from f and g states supports the accuracy of 
the recombination rates for this ion. That remarkable con- 
sistency can be explained by three factors concerning the 
observed high angular momentum states: those states are 
less dependent on atomic calculation assumptions because of 
their reduced non-hydrogenic effects, their populations are 
free from optical depth effects, and their fine-structure com- 
ponents are often blended in spectroscopic terms. In those 
blends, the effective transition rate of the term is a weighted 
average of several fine structure transition probabilities so 
that the total intensity of the blend is the sum of intensities 
of those components. Quantum calculations with different 
assumptions tend to agree on the transition probabilities for 
the more intense components of a term. The average tends 
to cancel out differences in the calculated ^4-values of the 
fine structure components while the summation gives more 
weight to the components with the larger intensities and 
lower observational errors. Measurement errors are also ex- 
pected to decrease with the higher S /N of a blend of several 
lines belonging to the same term. 

Recombinations produce more than half the intensity of 
most Cll lines, but we find an important fluorescence con- 
tribution to all the observed lines from s, p and d states. 
SBW did not find differences in the line profiles of those 
lines with profiles of recombination lines as should be indi- 
cated by their proposed criteria based on line profile differ- 
entiation. The fluorescence contribution is particularly im- 
portant in multiplets 3p 2 P°-4s 2 Si /2 AA 3918.97, 3920.68, 
4p 2 P°-6d 2 DAA 4637.63-4639.07, 3d 2 D-4p 2 P°AA 5889.28- 
5891.60, 4p 2 P°-5d 2 DAA 6257.18-6259.80, and 3s 2 Si /2 - 
3p 2 P°AA 6578.05, 6582.88, and it is mostly due to ab- 
sorptions from the ground state in the doublet 2p 2 P°- 
4s 2 SAA 636.994, 636.251, and the triplets 2p 2 P°- 
5d 2 DAA 560.239-560.439, 2p 2 P°-6d 2 DAA 543.258-543.445, 
followed by decays to term 4p 2 P°. This fluorescence exci- 
tation explains the enhanced intensities of those multiplets 
relative to the recombination lines, and the seemingly high 
C 2+ abundances derived from their intensities by SBW us- 
ing recombination rates alone. 

We added 13 lines from the SWBH survey with excited- 
core upper terms 2s2p( 3 P°)3p 2 P and 2s2p( 3 P°)3p 2 D and 



two lines from the 9f 2 F° ter m using addi t ional effec- 
tive recombination coefficients of iDavev et al.l (|2000h pub- 
lished on line at the CDS site. Term 2s2p( 3 P°)3p 2 P is 
mostly excited by fluorescence through absorptions in mul- 
tiplet 2p 2 P°-2s2p( 3 P°)3p 2 PAA 549.320-549.570. The cal- 
culated intensities for decays from that term agree with 
the observations within our adopted uncertainty intervals 
except for the line 2p 3 2 P° /2 -2s2p( 3 P°)3p 2 P 3/2 A 7508.89 
apparently undetected by SWBH, and the line 4p 2 P 3 / 2 - 
2s2p( 3 P°)3p 2 P 3/2 A 5121.83 underestimated by a factor of 
2.8. Lines from term 2s2p( 3 P°)3p 2 D on the other hand 
are underestimated by our calculations by 30 to 50 per 
cent. Lines from the state 2s2p( 3 P°)3p 2 D 3 / 2 have a 70 
per cent excitation by fluorescence while lines from the 
other state 2s2p( 3 P°)3p 2 D 5 / 2 have a 40 per cent exci- 
tation by fluorescence mainly in absorptions in multiplet 
2p 2 P°-2s2p( 3 P°)3p 2 DAA 530.275-530.474. The relatively 
intense line 4p 2 P° /2 -2s2p( 3 P°)3p 2 D 3/2 A 3835.72 is proba- 
bly blended with the Balmer H9 line. The two blended lines 
5d 2 D 5/2 -9f 2 F° 5/2 and 5d 2 D 5/2 -9f 2 F? /2 at A 7860.50 are 
the only f multiplet with a strong discrepancy with obser- 
vations with an observed intensity 2.8 times larger than the 
predicted intensity. The discrepancy may be due to uncer- 
tainties in the measurement because the observed S/N of 
28.7 in SWBH is lower than the other C n blends of f and 
g states and the intensity of the other blend from the 9f 2 F 
term is within the uncertainty interval of the observed value. 



5.4 Nil 

Fluorescence is the main excitation mechanism for most 
Nil lines in IC 418, but to a lesser extent than in the 
Orion nebula. The main p umping mechanisms in t his io n 
were discussed in detail bv lEscalante and Morissetl (2005). 
The 3 P ground core offers several pumping channels to 
the 3d 3 P° and 3d 3 D° terms through multiplets 2p 2 3 P- 
3d 3 P°AA 529.36-529.87, 2p 2 3 P-3d 3 D°AA 533.51-533.88 
and 2p 2 3 P-4s 3 P°AA 508.48-509.01. The only 3d term with 
significant recombination excitation is 3d 3 F°, but its bright- 
est components at AA 5001.13, 5001.13 and 5005.15 are 
blended with the [O iii]A 5006.84 nebular line. The singlet 
lines 3p ^i-Sd X D^A 4447.03 and 3p x D 2 -3d ^^6610.56 
listed by SBW are also excited mostly by fluorescence due to 
absorptions in the transitions 2p 3 Pi-3d X D 2 A 534.657 and 
2p 1 D 2 -3d X F 3 A 574.650 respectively. The first one is a spin- 
forbidden transition with the ground triplet term while the 
second demonstrates that a metastable state can also pro- 
duce fluorescence excitation. 

The disagreement of observed intensities with predicted 
intensities in lines excited by fluorescence is greater in this 
ion. There is a systematic underestimation of the calculated 
values of lines from 3d states of 30 to 70 per cent. The rest 
of the calculated intensities selected by SBW are within 40 
per cent of the observed values. 

The 4f states are populated by recombination. Direct 
recombinations and transitions from the 4f states produce 
between 30 and 60 per cent of the intensity of the lines from 
the 3d 3 D° and 3d 3 F° terms. EM showed that fluorescence 
cannot significantly populate f and higher angular momen- 
tum states in N + because quantum probabilities strongly 
favour downward transitions with decreasing angular mo- 
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Figure 2. Comparison of C II intensities observed by SBW with predicted intensities produced by recombination and fluorescence of 
lines from s states (black squares), p states (black triangles), d states (open red circles), and f and g states (filled blue circles). The 
dash at the bottom of the vertical lines marks the intensity produced by recombination excitation only. The broken horizontal line is the 
average value of / C al/^obs- Note that the observed intensities are normalized to /(H/9) = 10 4 . 
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Figure 3. Same as in Fig. [2] for Nil lines. 



mentum. That is a well known trend of transition probabil- 
ities inferred from atomic data compilations and theoreti- 
cal considerations (|Khriplovich fc Matvienkoll2007l ). Conse- 
quently d states have a greater probability of decaying to 
p states than to f states. Should fluorescence enhance an f 
state population by decays from a higher d state, lines from 



lower p states, like 3s 3 P2-3p 3 P 2 A 4630.54, would show at 
least a 70 per cent enhancement in their predicted inten- 
sity (EM). In our present calculations that line is already 
overpredicted by 20 per cent above its observed value with- 
out any plausible fluorescence contribution to f states. The 
enhanced intensity of the 3s 3 P2~3p 3 DsA 5679.56 line by 
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Figure 4. Same as in Fig. [2] for On lines. 



fluorescence in IC 418 puts its ratio with the recombination 
line 3d 3 F%-4i G[9/2] 5 A 4041.31 out of the T e and n e ranges 
in fig. 8 of FSL, and highlights the importance of determin- 
ing the atomic processes that excite a line before it is used 
as a diagnostic tool. 

We have added 15 lines from the SWBH survey 
that were omitted in the analysis of SBW. Some like 
3d 3 D§-4f F[5/2] 3 A 4241.76, and 3d 3 D§-4f F[7/2] 4 A 4241.78 
are blended, and 3d 3 Dg-4f F[7/2] 3 A 4237.05 is probably 
blended with the transition 3d 3 D?~4f F[5/2] 2 A 4236.91 
listed in the selection of SBW. There are four lines 
standing out from f states with predicted intensi- 
ties significantly below the observed values: 3d 3 F!j- 
4f G[7/2] 3 A 4058.16, 3d ^SHf F[5/2] 3 A 4176.16, 3d 3 D£- 
4f F[5/2] 2 A4246.71, 3d 3 Pg-4f D[3/2]i A4433.48, and 
4d ^3-5! G[9/2] 4 A 9794.05, which are discussed below. 

Because there is only one more line from a 4f D term 
with a secure identification (3d : P?-4f D[3/2] 2 A 4678.14) 
and an intensity 15 per cent above the predicted value, 
we cannot assert whether the effective recombination coef- 
ficient of the A 4433.48 line is underestimated. The other 
line produced by a 4f D term listed in SWBH, 3d 3 P 2 - 
4f D[5/2] 3 A 4432.74, may be blended with a [Fe II] line, but 
our calculations predict an intensity 30 per cent above the 
measured value. The effective recombination coefficients of 
FSL and EV at 10 4 K agree for the lines of multiplet 3d 3 P°- 
4f D [5/2] AA 4432.74-4442.02. 

Our calculations predict intensities 0.002 times and 
0.006 times the observed intensities of the lines A 4058. 16 
and A 4246.71 respectively. The low S/N of 10.3 and 9.2 and 
the displacement between the tabulated and observed wave- 
length of 10.3 km s -1 and 9.2 km s^ 1 respectively reported 
by SWBH for those lines suggest that they may be misiden- 
tifications. Other lines of the term 4f F[5/2], like 3d 1 D§- 



4f F[5/2] 3 A 4176.16, with more secure identifications also 
show important differences with the calculated values. 

The dispersion of calculated intensities above and below 
the observed values of transitions from 4f D and 4f F terms 
may imply uncertainties in the recombination coefficients or 
uncertainties in the measurements due to the low signal to 
noise ratio. There are no significant differences between the 
effective recombination coefficients of FSL and EV for those 
lines as shown in table [1] 

For the lines from 4f G terms the agreement between 
predicted and observed intensities is within 20 per cent for 
nearly all lines. The only exception is the line 3d 3 F 4 - 
4f G[7/2] 3 A 4058.16 mentioned above with an effective re- 
combination coefficient by FSL more than twice the value 
of EV. Even so, the value by FSL is not large enough to 
explain its intensity if it is indeed N n emission. 

Finally the line 4d 1 F 3 ~5i G [9/2] 4 A 9794.05 is the only 
line from a 5f G term above the detection limit, and is 
strongly underestimated by a factor of 7.6 in intensity, with 
a S/N of 9.0 and no alternative identifications. 

The line at A 4674.903 was measured but not iden- 
tified by SWBH. Our calculations predict the Nil spin- 
forbidden transition 3s 1 P?-3p 3 P A 4674.91 to be a 
m ost likely identification. The tran sition probability given 
by iFroese Fischer fe Tachievl (|2004h of 8.54 x 10 s over- 
estimates the observed value by 40 per cent. The other two 
members of this Nil multiplet, 3s ^J-Sp 3 P 2 A 4654.53 and 
3s ^l-Sp 3 Pi A 4667.21, were identified by SWBH as the 
Ol multiplet 3p 5 P-8d 5 D° A 4654.56 and the line [Fern] 
3d 6 5 D 3 -3d 6 3 F2 2 A 4667.010 respectively. Our calculations 
predict an intensity 0.09 times weaker than the observed 
value for the Ol multiplet (see section I5T5)) while the agree- 
ment with predicted intensities for both lines is much better 
if our proposed identifications as Nil lines are correct. 
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Table 5. Predicted and observed line intensities of Cn. /rec/^calc 
is the fraction of the calculated intensity due to recombination, 
^calc is the predicted intensity, / b s is the measured intensity with 
a 20 per cent minimum uncertainty. The ? marks measured in- 
tensities taken from SWBH, but not selected by SBW. 



Table 5 - continued 



Lower-Upper 

4p 2 P° /2 -3p' 2 D 5/2 t 
4p 2 P° /2 -3p' 2 D 3/2 
4p 2 P° /2 -3 P ' 2 D 3/2 
3p 2 P° /2 -4s 2 S 1/2 
3p 2 PS - A - 2c 



3/2 



-4s % /2 
3d 2 D 5/2 -4f 2 F° /2 
4f 2 F°-10g 2 G § 
4d 2 D-9f 2 F° 
4f 2 F°-9g 2 G 
4d 2 D-8f 2 F° 
4p 2 P° /2 -6d 2 D 3/2 
4p 2 P° /2 -6d 2 D 5/2 
4 P 2 P° /2 -6d 2 D 3/2 
2 P 2 2 P 3 /2-3 P a Pg /a 
4f 2 F°-8g 2 G 
2P 3 2 P? /2 -3 P ' 2 D 5/2 
2 D- 



2p 3 2 P? /2 -3p' 

,3 2 



2p 3 2 P°, -3p' 2 D 



3/2 

4p 2 P? /2 -3p 



J 3/2 
3/2 



3/2 



4p 2 P° /2 -3p' 2 P 3/2 
4p 2 P° /2 -3p' 2 P 1/2 



1/2 

4p 2 P° /2 -3p' " 
4f 2 F°-7g 2 G 



1/2 



3d 


2 D 3 /2- 


dp 


2po 
*3/2 


3d 


2 D 5/2 - 


dp 


2po 
r 3/2 


3d 


2 D 3 /2- 


dp 


2po 

1/2 


4d 


2 D-6f 


2 p. 


> 


dp 


2po 
r l/2" 


-5d 


2 D 3 /2 


4p 


2po 
r 3/2~ 


-5d 


2 D 5/2 


dp 


2po 
r 3/2" 


-5d 


2 F>3/2 



4f 2 F°-6g 2 G 
3s 2 S 1/2 -3p 2 P° /2 
3s 2 S 1/2 -3p 2 P° /2 
3p 2 P° /2 -3d 2 D 3/2 



1/2 
o 

3/2 



3p 2 P? /0 -3d 2 D 



5/2 



3p 2 P° /2 -3d 2 D 3/2 



J 9/2 



■ 3 /2 

5f 2 F° /2 -10g 2 G 7/2 
5f 2 F° /2 -10g 2 G 9 
2p 3 2 P° /2 -3 P ' 2 P 3/2 
2p 3 2 P° /2 -3p' 2 P 3/2 
2p 3 2 P° /2 -3 P ' 2 P 1/2 
2p 3 2 P° /2 -3p' 2 P 1/2 
5d 2 D 3/2 -9f 2 F° /2 
5d 2 D 5/2 -9f 2 F° /2 
5f 2 F° /2 -9g 2 G 7/2 



5f 2 F° /2 -9, 



% G, 



9/2 



5d 2 D 3/2 -8f 2 F° /2 



7/2 



5d 2 D 5/2 -8f 



A (A) 


-^rec /^calc 


-^calc 


383d. 7 


0.5d2 


2.d0 


3835.7 


0.226 


2 5d 


OOOO. { 




U.OU 


Qm n n 
oy ly .U 


n i no 


Q Q7 
O.O / 


QOOO 7 




i 7 7n 


A OK7 Q 


u.yoy 


c^ri on 


A OCIO O 


l.UUU 


U.ol 


4329.7 


0.954 


0.68 


449d.d 


d.000 


d.d8 


46d8.8 


0.942 


0.97 


4637.6 


0.492 


0.56 


4638.9 


0.807 


0.65 


4639. d 


0.492 


O.dd 


4744.8 


0.654 


0.09 


480d.9 


d.000 


d.82 


5032. d 


0.5d2 


2.48 


5035.9 


0.226 


2.97 


5040.7 


0.226 


0.59 


5d20.d 


0.077 


0.d9 


5d2d.8 


0.077 


0.96 


5d25.2 


0.043 


0.69 


5d27.0 


0.043 


0.34 


534d.5 


d.000 


3.dd 


5889.3 


0.506 


0.30 


5889.8 


0.506 


2.69 


589d.6 


0.406 


d.86 


6d45.9 


0.969 


2.22 


6257.2 


0.270 


0.67 


6259.6 


0.544 


0.64 


6259.8 


0.270 


0.d3 


6460. d 


d.000 


6.04 


6578. d 


0.654 


42.90 


6582.9 


0.527 


26.60 


723d. 3 


0.632 


26.60 


7236.4 


0.836 


43.30 


7237.2 


0.632 


5.30 


7508.9 


d.000 


0.d7 


7508.9 


d.000 


0.23 


7508.9 


0.077 


0.25 


75d9.5 


0.077 


d.23 


75d9.9 


0.043 


0.89 


7530.6 


0.043 


0.44 


7860.2 


0.943 


0.09 


7860.5* 


0.962 


0.d3 


8d39.6 


d.000 


0.24 


8d39.6 


d.000 


0.3d 


8867.6 


0.923 


O.dd 


8868. d 


0.956 


0.d6 



3.00 



d0.68 

20.52 

57.d2 

0.79 

0.70 

d.09 

d.07 

0.56 



d.83 
4.34 ? 
5.58 ? 



2.68 ? 
0.58 ? 
0.30 ? 
2.77 



d.8d 
2.53 
0.52 
0.84 

5.84 
53.74 

d6.92 
46.73 
4.89 



d.03 ? 
0.92 ? 
0.38 ? 

0.36 ? 



§ Terms with no subscript are blended transitions of several states 
of angular momentum J. 

t 3p' is core excited configuration 2s2p( 1 D°)3p. 
*A 7860.5 includes blend with 5d 2 D 5/2 -9f 2 F° /2 . 



Lower-Upper 


A (A) 


^rec/^calc 


-^calc 


^ohs 


5f 2 F° /2 -8g 2 G 9/2 
5f 2 F° /2 -8g 2 G 7/2 
3s' 2 P° /2 -3 P ' 2 D 5/2 
3s' 2 P° /2 -3 P ' 2 D 3/2 


922d.d 
922d.d 
9238.3 
925d.O 


d.000 
d.000 
0.5d2 
0.226 


0.45 
0.35 
0.d4 
0.d7 


0.24 ? 


(Irec/ -^calc) 


= 0.58d 


{-^calc 


/lobs) 


= 0.938 



5.5 On 

The lines in this ion are produced mostly by recombination, 
which strongly favours the quartets because of its higher 
multiplicity. Fluorescence is limited to excitations from the 
2p 3 4 S 3 / 2 state, and thus contributes less than half the inten- 
sity of lines from the 3p 4 P° term through absorptions in the 
2p 3 4 S°-3d 4 PAA 429.918-430.176 multiplet with some con- 
tribution from absorptions in the 2p 3 4 S°-3d 4 DAA 429.650- 
429.716 multiplet, followed by decays to the 3p terms. Most 
doublets in the ground core configuration are being excited 
by recombination with less than 20 per cent contribution 
of fluorescence from the 2p 3 2 D° metastable term. The ex- 

o 

3/2 



ception are the states of the 3d F term. Lines 3p D 



3d 



F 5/2 A 4699.22 and 3p 2 Dg /2 -3d 2 F 7/2 A 4705.35 have a 



term 
3d 



50 per cent contribution of fluorescence from the metastable 
2p 3 2 D°. The spin-forbidden transition 2p 3 4 S 3/2 - 
2 F 5/ / 2 A 429.560 contributes 10 per cent to the fluores- 
cence excitation of the A 4699.22 line. This line may be 
blended with the core excited transition 2p 2 ( 4 D)3p 2 Dg/ 2 - 
2p 2 ( 1 D)3d 2 F 7/2 A4699.011 for which we lack a published re- 
combination coefficient, and probably contributes to the in- 
tensity of the line, which appears somewhat underestimated 
in our calculations. 

The dielectronic reco mbination coefficients from 
iNussbaumer fc Storevl |l984) give a good fit to the ob- 
served intensities of lines of configurations with the excited 
core 2p 2 ( 1 D). We estimated the co e fficien t of some lines 
not listed by INussbaumer fc Storevl (jT984) by scaling the 
coefficients of other lines by the corresponding branching 
ratios as explained in section 15.11 Nevertheless many 
other O II lines with an excited core listed in SWBH lack 
effective recombination coefficients and were left out of our 
calculation. 

We added eight lines identified by SWBH to the selec- 
tion in SBW. Some may be blended with transitions from 
other ions, but O II probably contributes most of the inten- 
sity of the line. Line 3p 4 P° /2 -3d 4 D 3/2 A 4097.22 is probably 
blended with line 3d 4 F 7/2 -4f G[4] 9/2 A 4097.26 and we have 
added their calculated intensities in table [7] 

The largest discrepancy in our calculations occurs in the 
line 3p 4 P 3 / 2 -3d 4 D 1/ / 2 A 4110.79 with a calculated intensity 
a factor of 11 below the observed value. Lines with upper 
terms 3d 4 P and 3d 4 D also tend to be underestimated in 
our calculations, and the reason appears to be uncertainties 
in the transit ion probabilit ies. The A-values for those terms 
calculated bv lNaharl l|2009h are sig nificantly lower than t hose 
quoted in t he NI ST database from I Veres fc Wise! l|l996T ) and 
iBell et al.l (ll995T ). The NIST values give a much better - 
albeit insufficient- match to observations with predicted in- 
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Table 6. Predicted and observed line intensities of Nil. 



Lower-Upper 


A (A) 


^rec / A; ale 


-feale 


^obs 




dp 


3 c, 


Ac, 3po 


3593.6 


0.004 


1.79 






dp 




Ac 3do 
-4S r ^ 


3609.1 


0.003 


1.41 






Q,, 

dp 


3c, 


Ac 3 do 

-4s P Q 


3615.9 


0.007 


0.22 






dp 


3p 

n- 


Ac- 3 DO 

-4s F 2 


3829.8 


0.004 


3.36 


2.32 




dp 


3p„ 
r 2 


Ac- 3 DO 

-4s P 2 


3838.4 


0.004 


9.68 


7.34 


? 


■3 

op 


3p 

Po 


Ac- 3 DO 
-4S r ^ 


3842.2 


0.003 


2.87 


1.26 




■3 

OP 


3p 
rt 


Ac- 3 do 
-4S r ^ 


3847.4 


0.003 


2.08 






dp 


3p 


Ac. 3po 
-4s P Q 


3855.1 


0.007 


1.19 


0.98 




dp 


3p„ 


Ac- 3 DO 
-4S r ^ 


3856.1 


0.003 


3.47 


10.13 ? 


du 


3po 


Af \c\ /Ol 
-41 vj|y/ Zj4 


4026.0 


0.999 


0.33 


— 




du 


3 t?o 
b 2 


-41 VjL ' /^J3 


4035.0 


1.000 


0.54 


0.68 




du 


3po 


/if rn /ol 
-41 Cj[y/2J 5 


4041.3 


1.000 


0.92 


1.19 




OU 


3 1? o 


yip z^ 1 /Ol . 
-41 U[//2j4 


4043.5 


0.999 


0.41 


0.41 




ou 


3po 


/If pf7 /ol 

-41 L ' / 3 


4058.1 


1.000 


0.00 


0.43 


? 


OU 


3 17 o 


Af T?\ r 7 /Ol . 

-41 t[(/Z\4 


4095.9 


0.910 


0.09 


— 




OU 


1 Dio 
D 2 


Af r^rv /Ol 
-41 t 


4171.5 


0.944 


0.33 


— 




QA 
OU 


D 2 


Af DTc; /ol . 
-41 £ [0/2J 3 


4176.1 


0.927 


0.25 


0.64 




■3, J 

OU 


3n° 

U 3 


Af T\\K /Ol 

-41 D|5/2j3 


4179.6 


0.923 


0.09 


— 




OU 


3 T\o 


Af 17 re; /ol~ 
-41 £ [0/4 2 


4236.9* 


0.900 


0.58 


0.30 


? 


3rl 
OU 


3no 
U 2 


/tf dTk /o"|„ 
-41 F [D/Z\3 


4241.7* 


0.915 


1.04 


0.69 


? 


'3/1 
OU 


U 3 


-41 F [0/ Z\2 


4246.7 


0.874 


0.00 


0.31 


? 


■3^1 
OU 


3po 
F l 


Af DiTQ /Ol 

-41 D[d/2j2 


4427.2 


0.952 


0.15 


— 




■3^1 
OU 


f i 


/If DiTQ /ol 

-41 U[d/2Ji 


4427.9 


0.930 


0.08 


— 




'3/1 
OU 


3po 
F 2 


/if nix /oi„ 

-41 1J[0/^J3 


4432.7 


0.923 


0.47 


0.36 


? 


QA 
OU 


F o 


/If Di TQ /Ol .. 

-41 JJ[d/^Ji 


4433.4 


0.930 


0.11 


0.28 




'3/1 
OU 


3po 
F l 


/if "n /oi„ 

-41 D[D/ZJ2 


4442.0 


0.967 


0.10 


— 




■3 

op 


lp 
rt 


*3,-l 1 DlO 

-ou L>2 


4447.0 


0.177 


0.43 


0.39 




dp 


3n, 


3 DO 

-3d P 


4459.9 


0.070 


0.09 


0.25 




■3 

op 




■3^1 3 DO 

— ou r-^ 


4465.5 


0.034 


0.11 


— 




'3 r-« 

op 


3n„ 

U2 


oj 3po 
-ou 


4477.7 


0.034 


0.43 


0.61 




■3 

OP 


3n„ 
L>3 


Q^l 3 do 


4507.6 


0.117 


0.36 


0.51 




OU 




/if Tq /ol , 
-41 vjx|y/ zj4 


4530.3 


0.999 


0.36 


0.42 




Q^l 
OU 


1 DO 

F 3~ 


/if z^ 1 /oi . 
-41 C_r[/y2j4 


4552.5 


0.999 


0.28 


0.26 




o„ 
OS 


3po 
F l" 


Q n 3 D„ 

dp 1^2 


4601.5 


0.126 


2.92 


2.63 




'}., 
OS 


3 DO 

P (T 


Qv, 3d 
dp 


4607.2 


0.110 


2.79 


2.57 


? 


os 


3 DO 

F r 


Qr, 3D 

dp ir*i 


4613.9 


0.110 


1.93 


1.77 




OS 


3 "DO 

F l" 


o„ 3d 
OP PO 


4621.4 


0.086 


3.42 


2.64 




OS 


3 DO 

F 2~ 


Q„ 3 D„ 

dp Jr*2 


4630.5 


0.126 


9.54 


8.05 




os 


3 DO 

F 2~ 


Qr, 3D 
dp -fl 


4643.1 


0.110 


3.83 


3.44 




os 


1 do 

F r 


dp f2 


4654.5 


0.126 


0.24 


0.31 


? 


Ob 


1 DO 

F r 


o„ 3d 
dp Fi 


4667.2 


0.110 


0.20 


0.28 


? 


'),, 
OS 


1 DO 

F i" 


o„ 3d,, 
-dp J^o 


4674.9 


0.086 


0.30 


0.21 


? 


'3,1 
OU 


1 DO 

F i 


/If T~l \Q /Ol 

-41 U[d/^j2 


4678.1 


0.952 


0.12 


0.14 




'3,1 
OU 


1 DO 

F i 


/If Di rt; /ol ^ 
-41 D[D/^j2 


4694.6 


0.967 


0.17 






dp 


3n 


oj 3n° 
— OU 


4774.2 


0.112 


0.16 


0.22 




■3 

dp 


3n, 


■3^1 3 DiO 

-du 


4779.7 


0.048 


1.21 


1.79 




■3 

op 


3n 


Q^l 3 DlO 

-dU Ug 


4781.2 


0.409 


0.06 


0.13 




■3 

op 


3n„ 

L>2 


1A 3 DiO 

-du U9 


4788.1 


0.112 


1.20 


1.95 




3p 


3 D 2 


-3d 3 D° 


4793.6 


0.048 


0.37 






3p 


3 D 3 


-3d 3 D° 


4803.3 


0.409 


0.89 


2.61 




3p 


3 D 3 


-3d 3 D§ 


4810.3 


0.112 


0.23 


0.40 


? 


3p 


3 Sr- 


-3d 3 Pg 


4987.4 


0.070 


0.53 






3p 


3 Sr- 


-3d 3 PJ 


4994.4 


0.034 


3.29 


4.23 




3p 


3 Dr 


-3d 3 F° 


5001.1 


0.668 


1.19 






3p 


3 D 2 


-3d 3 F§ 


5001.5 


0.493 


2.75 






3s 


3po 
^0 


■3p 3 Si 


5002.7 


0.119 


0.58 






3p 


3 D 3 


-3d 3 F° 


5005.2 


0.929 


2.26 






3p 


3 Si- 


-3d 3 P° 


5007.3 


0.117 


2.53 






3s 


3po 


3p 3 Si 


5010.6 


0.119 


1.49 






3p 


3 D 2 


-3d 3 F° 


5016.4 


0.668 


0.20 






3p 


3D 3 


-3d 3 F§ 


5025.7 


0.493 


0.28 







*A 4236.9 includes blend with 3d 3 Dg-4f F[7/2] 7 A 4237.047 
A 4241.7 includes blend with 3d 3 D°-4f F[7/2] 9 A 4241.786 
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Table 6 - continued 



Lower-Upper 


A (A) 


/ -^calc 


I al 


-fobs 


3s 3 P°-3p 3 Si 


5045.1 


0.119 


2.31 


2.89 


3s ^"j^p 3 Si 


5073.6 


0.119 


0.17 


— 


3p 3 P -3d 3 P° 


5452.1 


0.034 


0.35 


0.42 


3p 3 Pi-3d 3 P° 


5454.2 


0.070 


0.22 


— 


3p 3 Pi-3d 3 P° 


5462.6 


0.034 


0.40 


0.45 


3p 3 Pi-3d 3 P° 


5478.1 


0.117 


0.14 


0.30 


3p 3 P 2 -3d 3 P° 


5480.1 


0.034 


0.51 


0.56 


3p 3 P 2 -3d 3 P° 


5495.7 


0.117 


0.70 


1.48 


3s 3 P°-3p 3 D 2 


5666.6 


0.208 


6.01 


4.14 


3s 3 Pg-3p 3 Di 


5676.0 


0.252 


2.37 


1.97 


3s 3 P^-3p 3 D 3 


5679.6 


0.263 


9.40 


6.74 


3s 3 Pf-3p 3 Di 


5686.2 


0.252 


1.55 


1.27 


3s 3 P^-3p 3 D 2 


5710.8 


0.208 


1.98 


1.36 


3s 3 Pg-3p 3 Di 


5730.7 


0.252 


0.11 


0.13 


3s ^J-Sp 3 D 2 


5747.3 


0.208 


0.54 


— 


3s ^J-Sp 3 Di 


5767.5 


0.252 


0.19 


— 


3p 3 P -3d 3 D° 


5927.8 


0.048 


1.24 


1.91 


3p 3 Pi-3d 3 D§ 


5931.8 


0.112 


1.65 


2.73 


3p 3 Pi-3d 3 D° 


5940.2 


0.048 


0.87 


1.39 


3p 3 P 2 -3d 3 D§ 


5941.7 


0.409 


1.26 


3.15 


3p 3 P 2 -3d 3 D° 


5952.4 


0.112 


0.49 


0.52 


3ro 3r\„ 
oo r 1 ^ ^p u 3 


U1U / .o 


n ^1 q 
u.ory 


21 




3d 3 F°-4p 3 Di 


6170.2 


0.423 


0.11 


0.09 


3d 3 Fg-4p 3 D 2 


6173.3 


0.453 


0.16 


0.23 


3s 3 P°-3p x Pi 


6379.6 


0.351 


0.08 




3s ^-Sp 1 P 1 


6482.0 


0.351 


0.39 




3p 1 D 2 -3d J F° 


6610.6 


0.226 


0.17 


0.28 


3d 3 P§-4p 3 Si 


6810.0 


0.223 


0.09 


0.21 


3d 3 PJ-4p 3 Si 


6834.1 


0.223 


0.06 


0.16 


4p 3 D 3 -5s 3 P§ 


8676.1 


0.031 


0.08 


0.12 


4p 3 D 2 -5s 3 P° 


8699.0 


0.017 


0.07 


0.08 ? 


4d iFv-Sf G[9/2] 4 


9794.2 


0.999 


0.02 


0.16 ? 


(Irec/ -^calc) 


= 0.348 


{-^calc 


/lobs) 


= 0.863 



tensities 0.09 to 0.4 times the measured value, far below the 
assumed observational uncertainty. Some comparisons with 
the observed and calculated line intensities with different 
/1-values are given in table 

5.6 The velocity field 

The pumping transitions in fluorescence excitation will be 
Doppler shifted in an expanding nebula, and absorb the 
stellar flux at longer wavelengths. Our high-dispersion SED 
gives the chance to notice the effect of the velocity field of 
the gas on the intensities of some lines mostly excited by flu- 
orescence. Several authors have used diff erent velocity fields 
to model forbidden-line profiles in PN (|Neiner et al.lll996l : 
iGesicki et al.| [l996; Zhang 200j|). We have assumed two ac- 
celerating velocity fields, V oc R, and V oc R 4 with a max- 
imum velocity of 40kms _1 at the outer edge of the nebula, 
which are similar to fields considered in the literature. Ta- 
bles [5] to [7] were computed with the V oc R 4 field. Although 
PN shells are expected to expand, we also considered a static 
field with zero velocity for comparison purposes. 

Turbulence is another component of the velocity field 
that must be considered. Turbulence of 14kms _1 was 
needed to model broad line profiles in a PN with a Wolf- 
Rayet centra l star, but was not needed fo r O-st ar PNe by 
iNeiner et all l|l99rj) . iMorisset fc Stasinskal l|2006j) . however, 
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Table 7. Predicted and observed line intensities of On. 



Table 7 - continued 



Lower-Upper 



3p 4 D° /2 -3d 4 D 3/2 
3p 4 D° /2 -3d 4 D 5/2 
3p 4 D° /2 -3d 4 D 3/2 
3p 4 D° /2 -3d 4 D 7/2 
3p 4 D° /2 -3d 4 D 5/2 
3p 4 D° /2 -3d 4 D 3/2 
3p 4 D° /2 -3d 4 D 7/2 
3p 4 D° /2 -3d 4 D 5/2 
3p 4 D° /2 -3d 4 P 5/2 
3d 4 F 7/2 -4f F[3]° /2 
3d 4 F 9/2 -4f F[4]° /2 
3p 4 D° /2 -3d 4 F 3/2 
3p 4 D° /2 -3d 4 F 5/2 
3d 4 F 7/2 -4f G[5]° /2 
3p 4 D° /2 -3d 4 F 7/2 
3p 4 D° /2 -3d 4 F 9/2 
3p 4 D° /2 -3d 4 F 3/2 
3d 4 F 5/2 -4f G[4]° /2 
3p 4 D° /2 -3d 4 F 5/2 
3d 4 F 3/2 -4f G[3]° /2 
3d 4 F 9/2 -4f G[5]° 1/2 
3p 4 D° /2 -3d 4 F 7/2 
5/2 -4f G[3]° /2 

o 

9/2 
o 

5/2 



3d 4 F 

3d 4 F 7/2 -4f G[4]° 
3d 4 F 3/2 -4f D[3] 
3p 4 P° /2 -3d 4 D 5/2 
3p 4 P° /2 -3d 4 D 3/2 



7/2 



3d 4 F 5/2 -4f D[3]2 



3p 4 P° /2 -3d 4 D 1/2 
3p 4 P° /2 -3d 4 D 7/2 
3p 4 P° /2 -3d 4 D 5/2 
3p 4 P° /2 -3d 4 P 1/2 
3p 4 P° /2 -3d 4 P 1/2 
3p 4 P? /2 -3d 4 P 3/2 
3p 4 P° /2 -3d 4 P 3/2 
3p 4 P° /2 -3d 4 P 5/2 
3p 4 P° /2 -3d 4 P 3/2 
3p 4 P° /2 -3d 4 P 5/2 

-3d' 2 G 7/2 f 



3p' 2 F° /2 
3p' 2 F° /2 -3d' 2 G 9/2 
3d 4 D 7/2 -4f F[4]° /2 
3d 4 D 3/2 -4f F[3]° /2 
3d 4 D 5/2 -4f F[3]° /2 
3d 4 D 5/2 -4f F[3]° /2 
3d 4 D 1/2 -4f F[2]° /2 
7/2 -4f F[3]° /2 



3d 4 D 

3d 4 P 5/2 -4f D[2]° /2 
3d 4 D 3/2 -4f F[2]° /2 
3d 4 D 3/2 -4f F[2]° 



3d 2 F 



5/2 



-4f F[3] 



3d 4 P 5/2 -4f G[3]° 



3/2 

□ 

7/2 



7/2 



A (A) 


-^rec/-^calc 


^calc 


^obs 


3842.8 


0.605 


0.15 




3850.7 


0.686 


0.16 




3851.1 


0.605 


0.30 




3863.6 


0.825 


0.10 




3864.5 


0.686 


0.54 




3864.8 


0.605 


0.20 




3882.3 


0.825 


0.86 


0.63 


3883.1 


0.686 


0.10 


0.27 


3907.6 


0.711 


0.09 


0.35 


4048.2 


1.000 


0.10 




4062.8 


1.000 


0.19 




4069.6 


0.921 


1.46 


2.04 


4069.9 


0.924 


2.32 


1.99 


4071.2 


1.000 


0.26 




4072.3 


0.952 


3.41 


3.27 


4075.9 


0.956 


4.90 


4.42 


4078.8 


0.921 


0.53 


0.56 


4083.8 


1.000 


0.44 


0.49 


4085.2 


0.924 


0.65 


0.76 


4087.1 


1.000 


0.42 


0.45 


4089.2 


1.000 


1.54 


1.14 


4093.0 


0.952 


0.46 


0.32 


4095.6 


1.000 


0.30 


0.42 


4097.2* 


0.795 


1.96 


1.15 


4098.2 


1.000 


0.25 




4104.8 


0.686 


1.97 


1.60 


4105.1 


0.605 


0.76 


0.46 


4107.0 


1.000 


0.18 




4110.9 


0.868 


0.07 


0.75 


4119.2 


0.825 


2.23 


2.19 


4120.2 


0.686 


0.15 


0.68 


4121.3 


0.465 


0.41 


1.66 


4129.3 


0.465 


0.64 




4132.7 


0.655 


0.24 


0.83 


4140.8 


0.655 


0.42 




4153.4 


0.711 


0.70 


1.84 


4156.4 


0.655 


0.70 


1.59 


4169.2 


0.711 


0.96 




4185.5 


0.997 


0.79 


0.81 


4189.9 


0.997 


0.82 


1.24 


4275.5 


1.000 


0.84 


0.65 


A 07^ O 

4z / o.y 


l.UUU 


U. ID 




4276.3 


1.000 


0.12 




4276.8 


1.000 


0.31 


0.79 


4277.4 


1.000 


0.21 


0.33 


4277.8 


1.000 


0.15 




4281.3 


1.000 


0.08 


0.19 


4282.9 


1.000 


0.24 


0.26 


4283.6 


1.000 


0.15 




4285.6 


1.000 


0.29 


0.31 


4291.2 


1.000 


0.24 





*A 4097.2 includes emission of 3p 4 P° /2 -3d 4 D 3/2 

f 3p' and 3d' are core excited configurations 2p 2 ( 1 D)3p and 
2p 2 ( x D)3d respectively. 



Lower-Upper 


A (A) 


-Tree 


/-^calc 


-^calc 


-^obs 


3d 2p 5/2- 


AC TTirolo 
- 4i F W 5 /2 


4292.2 


1 


.000 





.14 






3d 4 P 3 / 2 


-4f D[2]° /2 


4294.8 


1 


.000 





.36 





.59 


3d U 7 /2 


-4f G[5] g/2 


4303.5 


1 


.000 





.09 






3d 4 P 5 /2- 


A f "P\ To! 

-4f D[3]° /2 


4303.8 


1 


.000 





.63 





.66 


3d 4 Pi/ 2 - 


A C T\ To! O 

-4f D[2] 3/2 


4307.3 


1 


.000 





.16 


1. 


.19 


3d 2 F 7/2 - 


AC T7> f A 1 O 

" 4f F M7/2 


4312.2 


1 


.000 





.10 






3d 2 F 7/2 - 


AC T71 f A I O 

"4f F[4] 9/2 


4313.5 


1 


.000 





.19 






o„ 4d 

3s Pi/2- 


3p 4 P° /2 


4317.0 





.605 


1 


.58 


1 


.59 


3d 4 P 3 / 2 - 


-4f D[3]° /2 


4317.7 


1 


.000 





.11 






3s P3/2- 


■3P 4P 5/2 


4319.6 





.747 


1 


.25 





.83 ? 


3s 4 Pi/ 2 - 


* ^1/2 


4325.8 





.522 





.31 






3d L»5/2" 


-4f G[4] 7/2 


4331.2 


1 


.000 










3d 4 D 7 / 2 


A £ r A 1 O 

-4f G[4]° /2 


4332.6 


1 


.000 





.15 






3s P 3 / 2 - 


3 P 4p S/2 


4336.7 





.605 





.67 


0. 


.37 


3d 2 F 5/2 - 


A £ f~\ \ A I O 

-4f G[4] 7/2 


4340.3 


1 


.000 





.31 






3d 2 F 7/2 - 


a £ s~\ ftri 

-4f G[5J (J/2 


4342.1 


1 


.000 





.85 






3d 4 D 5/2 


-4f G[3] 7/2 


4344.4 


1 


.000 





.18 






3s P 3 / 2 - 


-3p 4 P° /2 


4345.6 





.522 


1 


.92 


1. 


.92 


3s D 3/2 


-3p 2 D° /2 f 


4347.5 





.649 





.57 






3s P5/2- 


3 P 4P 5/2 


4349.5 





.747 


3 


.51 


3 


.30 


3s z D 5/ / 2 


-3P 2pi 5/2 


4351.2 





.571 





.69 


0. 


.80 


3d 2 F 5/2 - 


-4f G[3]° /2 


4353.5 


1 


.000 





.16 






o„ 4 T) 

3s P 5 / 2 - 


3P 4P 3/2 


4366.8 





.605 


1 


.67 


1. 


.36 


3d 2 F 7/2 - 


A £ f~\ \ A I O 

-4f G[4] g/2 


4371.7 


1 


.000 





.15 






Q„ 2r> 

3s P 3 / 2 - 


-3p 2 D° /2 


4414.9 





.812 


5 


.26 


3 


.03 


3s P1/2- 


3P D 3 /2 


4416.8 





.815 


2 


.95 


1 


.9 1 


2td 

3s P3/2- 


3P 2pi 3/2 


4452.2 





.815 





.54 





.88 ? 


o , 2t> 

3d ^P 3 / 2 - 


-4f D[2]° /2 


4466.4 


1 


.000 





.14 





.29 


3d 2 P 3 / 2 


A £ f~i fol O 

-4f G[3]° /2 


4477.9 


1 


.000 





.13 






3d 2 Pi /2 


-4f D[2]° /2 


4489.4 


1 


.000 





.10 





.07 


3d 2 P 3/2 - 


-4f D[3]° /2 


4491.2 


1 


.000 





.21 






3s U 5 / 2 


_3 P F 7/2 


4590.9 





.729 






1 

1 


.to 


3s D 3 / 2 


- 3 P' 2F 5/2 


4596.2 





.774 





.53 





.94 


3d D 3/2 


A £ T7> Tol O 

-4f F[3] 5/2 


4602.1 


1 


.000 





.27 


0. 


.22 


3d ^D 5/2 


a £ TTi r a 1 
-4f F[4]° /2 


4609.3 


1 


.000 





.66 





.60 


3d z D 3/2 


-4f F[2]° /2 


4610.2 


1 


.000 





.21 






3s 4 Pi/ 2 - 


-3p 4 D° /2 


4638.8 





.731 


2 


.07 


1 


.96 


3s 4 P 3 /2- 


"3P 4D 5/2 


4641.6 





.794 


1 


.68 


4. 


.22 


3s 4 P 5 / 2 - 


-3p 4 D° /2 


4649.2 





.795 


8. 


.35 


6 


.70 


3s 4 Pi /2 - 


■3p 4 D°. 


4650.9 





.688 


2 


.05 


2 


.17 


3s 4 P 3/2 - 


3P 4 D^ /2 


4661.6 





.731 


2 


.30 


2 


.31 


3s 4 P 3/2 - 


-3p 4 D° /2 


4673.8 





.688 





.35 





.41 


3s 4 P 5/2 - 


3P 4 D5/ 2 


4676.2 





.794 


1 


.59 


1 


.58 


3s 4 P 5/2 - 


3P 4 D 3/2 


4696.4 





.731 





.17 





.27 


3p 2 D° /2 


-3d 2 F 5/2 


4699.3 





.441 





.70 


1 


.33 ? 


3 P 2 °5/2 


-3d 2 F 7/2 


4705.2 





.472 


1 


.36 


1 


.83 


3p 4 S° /2 - 


-3d 4 D 5/2 


4856.4 





.686 





.11 






3 P 4S 3/2" 


-3d 4 P 1/2 


4890.8 





.465 





.41 


1 


.37 


3 P 4S 3/2" 


-3d 4 P 3/2 


4906.9 





.655 





.71 





.42 


3P ^3/2- 


-3d 4 P 5/2 


4924.7 





.711 





.41 





.89 



<W4alc} = °- 83 8 



<4alc/4b S > = 0.898 



f 3p' and 3d' are core excited terms 2p 2 ( 1 D)3p and 2p 2 ( 1 D)3d 
respectively. 
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Table 8. Comparison of pred icted and obs erved line intensities 
of On with A-values from m lNaharl d2009f> and (2) NIST. 



in a 3 km s 1 turbulent field with respect to a field without 
turbulence. 



Lower— Upper 



A (A) 



A(10 6 
(1) 



- 1 ) 
(2) 



-^calc/-^obs 

(1) (2) 



3p 


U 5/2 


3d 


4 P 5 /2 


3p 


4po 
r 3/2 


3d 


4 D 1/2 


3p 


4po 
1/2" 


3d 


*l/2 


3p 


4po 
r 3/2 


3d 


^5/2 



8.64 0.059 



0.254 



3 3 /2 -3d 4 D 1/2 4110.9 0.75 0.0321 80.5 0.00004 0.090 

4121.3 1.66 6.82 56.0 0.038 0.245 

4153.4 1.84 1.09 72.8 0.007 0.382 



have shown that broadened line profiles in PNe can be ex- 
plained with a non-spherical geometry with negligible tur- 
bulence. In our case turbulence increases the pumping prob- 
ability because it increases the width of the Voigt profile of 
a pumping transition (see equations [3] and [4] of EM). We 
considered a moderate turbulent field of 3kms _1 . Table [9] 
shows some of the lines that exhibit the largest variations in 
intensity with different velocity fields. 

The effect of a non-zero expansion velocity is stronger 
in lines that are pumped by a single transition from the 
ground or a metastable state like the Nil lines from the 
states 3d 3 Pg, 3d 1 D§ and 3d x Ft The intensities of those 
lines generally change by amounts larger than the assumed 
observational uncertainty between 20 and 80 per cent in the 
velocity field that we have adopted with respect to a static 
gas, and tend to match the observations better than with a 
static field. The Nil 3d ^^4447.0 line increases by 

a factor of 2 with respect to a static field because the line is 
mostly pumped by the transitions 2p 3 Pi-3d ^^A 534.657 
and 2p 2 *D 2 -3d 1 T)'^X 582.156, which lie next to peaks of the 
SED at longer wavelengths (see Fig. [TJ. 

Lines that are pumped by more than one transition gen- 
erally do not show large changes because variations in the 
intensity of the stellar flux at the wavelengths in the rest 
frame of the absorber tend to cancel out for the different 
transitions in a multiplet. One exception are lines from the 
Nil 3d 3 P° state with two main pumping transitions from 
the ground term at 2p 2 3 P -3d 3 P° A 529.355 and 2p 2 3 Pi- 
3d 3 P° A 529.491, which receive an increased stellar flux si- 
multaneously at a redshifted wavelength. 

The sensitivity of lines affected by fluorescence on the 
velocity gradient depends on the ion spatial distribution 
within the nebula. The variation of intensities between the 
static and accelerating fields with no turbulence is compara- 
ble or slightly larger than our assumed uncertainty of 20 per 
cent for most of the N II and O II lines in table Since the 
+ ion is more extended towards the inner parts of the neb- 
ula, O II lines are more sensitive to the field gradient, with 
variations in intensity larger than the assumed uncertainty 
in the V oc R field. The C + and N + ions are more con- 
centrated towards the outer parts. Their fluorescence lines 
are nearly insensitive to the choice of velocity field with the 
exception of the three Nil lines 3p 3 Di-3d 3 Pq A 4459.94, 
3p ^i-Sd X D^A 4447.03 and 3p x D 2 -3d ^^6610.56 from 
the states mentioned above. 

Turbulence is the most important factor in the intensity 
of the C II lines because of the high C + column density and 
optical depth of its pumping transitions. Most C II lines af- 
fected by fluorescence increase their intensity by 30 per cent 



5.7 Nl 

Fluorescence is important in atomic nitrogen because sev- 
eral resonant lines are just below the Lyman limit. The ob- 
served N I dipole-allowed lines are quartets in the 3s-3p ar- 
ray mainly produced by absorptions in multiplets 2p 3 4 S°- 
3d 4 P AA 953.415-953.970 and 2p 3 4 S°-4s 4 PAA 963.041- 
964.626 followed by decays to the 3p terms. SBW noted 
that the similarity of the line profiles of quartets and the 
N I forbidden lines points to fluorescence as the main excita- 
tion mechanism of the Nl quartets. A Skms" 1 turbulence 
increases the intensity of the Nl quartets by a factor of 3 
because of the large optical depth of the pumping transi- 
tions, and brings them within the uncertainty interval of 
the observed intensities. 

Fluorescence contributes 50 per cent of the forbidden 
4 S3 /0 - 2 D° /9 A 5200.26 and most of the forbidden 

,3 4qo ' 



line 2p o 3/2 - ^ 5/2 , 

line 2p 3 4 S3/ 2 - 2 D3/ 2 A 5197.90 due to absorptions from the 



metastable term 2p D° and to the unusually intense spin- 
forbidden transition 2p 3 4 Sg /2 -3d 2 F 5/2 A 954.104. The state 
3d 2 F 3/ / 2 preferentially decays to state 2p 3 2 H^^ 2 besides 
other doublets producing half of the total fluorescence con- 
tribution to the forbidden lines. If fluorescence is turned 
off, however, electron collisions alone can substitute the 
role of fluorescence in the excitation of the forbidden lines 
A 5200.26 and A 5197.90. The reason is that the population 
of the ground state increases with the lack of absorptions 
to excited states, and thus the rate of collisional excita- 
tions also increases. The other forbidden line, 2p 3 4 Sgy 2 - 
2 PJ/2 i3 / 2 AA 3466.50, 3466.54, is more heavily influenced by 
fluorescence, but it is beyond the red end of the spectrum 
of SWBH. 

The fluorescence contribution to the forbidden lines de- 
pends on the A-value of the spin-forbidden A 954. 104 men- 
tioned above. Intermediate calculations by Irlibbert et al.l 
(Il99ll) gave a value of 7.609 x 10 5 s 1 , but the calculation 
bvlFroese Fischer fc Tachievl |2004h gave 1.954 x 10 7 s _1 and 
iGoldbach et all (|l992l ) measured 3.30 x 10 7 s 1 , which is the 
value listed in the NIST database and adopted by us. 

The measured intensities of the forbidden lines show 
the largest discrepancy with our predicted values for Nl. 
The agreement with the fluorescence lines improves if the 
density and size of the PDR are increased, but the forbid- 
den lines remain underestimated for reasonable values of 
the PDR parameters in equation (J?}. As shown by EM in 
Orion, the efficiency of fluorescence excitation grows in a 
manner analogous of the curve of growth. An increasing col- 
umn density eventually saturates the pumping transitions 
and the fluorescence excitation grows logarithmically in the 
line wings. 

We only added the two lines 3s 4 Pi /2 -3p 4 P3 /2 A 8188.01 
and 3s 4 P 3/2 -3p 4 P3 /2 A 8210.72 from SWBH to the selection 
in SBW. The multiplet 3p 4 D°-3d 4 DAA 9776.90-9872.15 
has three intense lines detected by SWBH. Its other de- 
tectable components at longer wavelengths are beyond the 
observed spectral range. 
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Table 9. Variation of normalized intensity of selected lines with velocity field. (1) V = and no turbulence, (2) V oc R and no turbulence, 
(3) V oc R 4 and no turbulence, (4) V oc R and 3 km s^turbulence, (5) V oc R 4 and 3 km s _1 turbulence. 



Lower-Upper 


A (A) 


-^calc/^obs 


-^obs 






(1) (2) (3) (4) 


(5) 


Cn 



3p 2 P° 1/2 -4s 2 S 1/2 
3p 2 P° /2 -4s 2 S 1/2 
4p 2 P° /2 -3p' 2 P 1/2 
2p 3 2 P° /2 -3p' 2 P 1/2 

2 P 3 2p 3/ 2 - 3 P' 2p V2 
4p 2 P° /2 -3p' 2 P 3/2 

2p 3 2 P° /2 -3 P ' 2 P 3/2 

4p 2 P° /2 -3p' 2 D 5/2 

2 P 3 2P 3/2- 3 P' 2 °5/2 

4p 2 P° /2 -5d 2 D 3/2 
4p 2 P° /2 -5d 2 D 5/2 



3919.0 


0.64 


0.64 


0.63 


0.85 


0.83 


10.6< l 


i 


3920.7 


0.66 


0.66 


0.65 


0.88 


0.86 


20.52 


5127.0 


0.87 


0.88 


0.87 


1.16 


1.14 


0.30 


? 


7519.9 


0.73 


0.91 


0.73 


0.98 


0.97 


0.92 


? 


7530.6 


0.88 


0.90 


0.89 


1.18 


1.16 


0.38 


? 


5121.8 


0.28 


0.27 


0.28 


0.34 


0.35 


2.68 


? 


7519.5 


0.93 


0.74 


0.73 


1.16 


1.20 


1.03 


? 


3831.7 


0.53 


0.78 


0.62 


0.94 


0.70 


3.00 


? 


5032.1 


0.43 


0.64 


0.50 


0.76 


0.57 


4.34 


? 


6257.2 


1.14 


0.94 


1.10 


1.08 


1.29 


0.52 




6259.6 


0.70 


0.60 


0.68 


0.65 


0.76 


0.84 





Nil 



3p 


3 D r 


3d 


3po 

r o 


4459.9 


0.43 


0.23 


0.32 


0.25 


0.35 


0.25 


3p 


3 D 2 


-3d 


3po 
*1 


4477.7 


0.48 


0.61 


0.62 


0.68 


0.70 


0.61 


3p 


3 Si- 


3d 


3po 
r l 


4994.4 


0.54 


0.68 


0.68 


0.76 


0.78 


4.23 


3p 


3 P 2 - 


-3d 


3po 
r 1 


5495.7 


0.46 


0.34 


0.41 


0.38 


0.48 


1.48 


3p 


3 Dr 


-3d 


3 D° 


4779.7 


0.66 


0.47 


0.62 


0.51 


0.68 


1.79 


3p 


3 Po- 


3d 


3 D° 


5927.8 


0.63 


0.45 


0.60 


0.50 


0.65 


1.91 


3p 


x Pr 


3d 


1 r>° 


4447.0 


0.68 


1.40 


1.12 


1.37 


1.11 


0.39 


3p 


X D 2 


-3d 


lpo 
r 3 


6610.6 


0.50 


0.78 


0.62 


0.77 


0.61 


0.28 



On 



3p 


4po 

A 1/2 ' 


3d 


4 Pl/2 


3p 


4qo 
°3/2 


-3d 


4 Pl/2 


3p 


4 r>° 


-3d 


^5/2 


3p 


4po 
^3/2" 


3d 


4 P 5 /2 


3p 


°3/2 


-3d 


4 P 5 /2 


3p 


4po 

r 3/2' 


3d 


4D 3/2 


3p 


2 D o . 
^3/2 


-3d 


2 F 5/2 


3p 


2 D° 
^5/2 


-3d 


2 F 7/2 



4121.3 


0.23 


0.35 


0.23 


0.38 


0.24 


1.66 


4890.8 


0.28 


0.43 


0.28 


0.46 


0.30 


1.37 


3907.6 


0.25 


0.35 


0.25 


0.34 


0.25 


0.35 


4153.4 


0.38 


0.52 


0.38 


0.52 


0.38 


1.84 


4924.7 


0.46 


0.63 


0.46 


0.62 


0.46 


0.89 


4105.1 


1.47 


1.40 


1.44 


1.59 


1.65 


0.46 


4699.3 


0.57 


0.60 


0.54 


0.60 


0.52 


1.33 ? 


4705.2 


0.64 


0.97 


0.77 


0.96 


0.74 


1.83 



5.8 Ol 

The triplet terms are excited by fluorescence in the PDR 
of IC 418 through a large number of resonant transitions 
of comparable strength from the grou nd term 2p 2 3 P, a sit- 
uation also encountered in Orion by iLucvl |2002r i. Triplet 
1 line profiles shown by SBW are similar to those of the 
1 forbidden lines in support of the fluorescence excitation 
mechanism. The quintets on the other hand are produced in 
the ionized region by recombination of + with no contribu- 
tion from fluorescence because they are not connected with 
strong transitions from the ground term and recombination 
preferentially populates higher multiplicity states. 

SWBH identified several multiplets between terms 
3p 3 P and 3p 5 P and s states with principal quantum num- 
ber n — 5 to 8 and d states with principal quantum number 
n = 4 to 7. SBW does not list them although their iden- 
tifications seem cert ain and only thei r effec tive recombina- 
tions are lacking in iPecmignot et all |l991). The effective 
rec ombination coefficients fo r those terms were calculated 
bv lEscalante fc Victorl l|l992T ). Our calculations tend to un- 



derestimate systematically the quintet and some triplet in- 
tensities by a factor of 1.3, which is slightly higher than the 
assumed uncertainty an d suggests that the e ffectiv e recom- 
bination coefficients of lEscalante fc Victorl (|l992T ) may be 
underestimated. The e ffective recombination coefficients of 
IPecm ignot et all (Il99ll ) are higher by about that same fac- 
tor, and thus give a closer agreement with the observations. 

The forbidden lines of the ground configuration 
2p 4 3 P 2 -2p 4 ^aA 6300.30, 2p 4 3 Pi-2p 4 ^aA 6363.78 and 
2p 4 1 D 2 -2p 4 x Si A 5577.34 are excited by electron collisions 
with negligible contribution from fluorescence in the warm 
interface of the PDR and ionized region. Our predicted in- 
tensities for those lines are underestimated by the same 
amount as the recombination lines, but they are more sen- 
sitive to the density than the Nl optical forbidden lines. If 
the ri3 parameter in equation (|7J) is raised to 35000, our pre- 
dicted intensities of the three forbidden lines perfectly match 
the observed values. This value of the density is already too 
high to explain t he FIR [C n] and [0 1] line observations of 
iLiu et all l|2001al ). It is also possible to improve the agree- 
ment of the calculated fluorescence and recombination lines 
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Figure 5. Same as in Fig. [5] for Nl lines. Collisionally excited forbidden lines are also included. 
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Figure 6. Same as in Fig. [5] for Ol lines. 



with the measurements if the PDR were larger, but there is 
no evidence of a dense and extended neutral medium around 
IC 418. Such a medium would increase the FIR line fluxes 
beyond their observed uncertainty. 

The largest discrepancy with the observations occurs 
in multiplet 3p 3 P ,i,2-5s 3 S?A 7254.15-53 with an ob- 
served intensity 1.8 times the predicted one. The possibility 
of an overestimated observed intensity due to a low S/N 



jRola fc Pelalll994) can be ruled out because the S/N mea- 
sured by SWBH for the A 7254 multiplet is 15.0. There is 
only one O I line with S /N < 7 in our calculations. Most of 
the O I lines or blends are relatively intense with S/N > 100. 
Lines from other states in the Rydberg series, 6s 3 S?, and 
8s 3 Si agree better with the calculated values. 

SWBH also identified the feature at A 4654.475 as the 
Ol multiplet 3p 5 P 2 -8d 5 D? >2 ,3 < 4654.56, which is strongly 
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Table 10. Predicted and observed line intensities of Nl. 



Table 11. Predicted and observed line intensities of O I. 
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OS 


4 Pi/2- 


-3p 


4co 
5 3/2 


'3.-. 
OS 


4 P 3 /2- 
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^1/2" 


-3p 
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3s 
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3s 
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-3p 


U 3/2 


3s 
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-3p 


4 D° 
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3s 


4 P 3 /2- 


-3p 


4 n° 

U l/2 


3s 


4 P 5 /2- 


-3p 


4 D° 
^3/2 



A (A) 



=/4. 



^calc 



3 P 4 D° /2 - 



J 3/2 



3p 4 D° /2 -3d 4 D 5/2 
-3d 4 D 1/2 



3p 4 D ?/2 



3p 4 D° /2 -3d 4 D 3/2 
3p 4 D° /2 -3d 4 D 1/2 



3p 4 D° /2 -3d 4 D 



5/2 



3p 


4t->o 

U 5/2 


-3d 


4 D 3 /2 


3p 


4 n° 

U 7/2 


-3d 


4 D 5/2 


3p 


4po 
r l/2 


-3d 


4 D 3/2 


3p 


4po 
*3/2 


-3d 


4 D 5/2 


3p 


4po 

A 1/2 


-3d 


4 D 1/2 


3p 


4po 
^3/2 


-3d 


4 D 3 /2 


3p 


4po 
r 3/2 


-3d 


4 D 1/2 


3p 


4po 
r 5/2~ 


-3d 


4 D 5/2 


3p 


4po 
r 5/2~ 


-3d 


4 D 3/2 



5197.9 
5200.3 
7423.6 
7442.3 
7468.3 
8184.9 
8188.0 
8200.4 
8210.7 
8216.3 
8223.1 
8242.4 
8680.3 
8683.4 
8686.1 
8703.2 
8711.7 
8718.8 
8728.9 
8747.4 
9776.9 
9786.8 
9788.3 
9798.6 
9810.0 
9822.7 
9834.6 
9872.1 
10500.3 
10507.0 
10513.4 
10520.6 
10533.8 
10549.6 
10563.3 



0.000 
0.000 
0.007 
0.007 
0.007 
0.017 
0.007 
0.005 
0.007 
0.017 
0.005 
0.007 
0.041 
0.022 
0.015 
0.013 
0.015 
0.022 
0.013 
0.015 
0.003 
0.005 
0.002 
0.003 
0.002 
0.005 
0.003 
0.005 
0.003 
0.005 
0.002 
0.003 
0.002 
0.005 
0.003 



3.58 
2.21 
1.62 
3.36 
5.21 
1.81 
4.36 
1.20 
1.66 
4.69 
6.37 
4.64 
3.40 
3.50 
2.12 
2.36 
2.48 
1.31 
0.42 
0.20 
0.19 
0.17 
0.42 
0.44 
0.74 
0.76 
0.72 
0.45 
0.92 
1.79 
2.37 
2.29 
0.95 
1.64 
0.49 



20.11 

11.73 

1.56 

3.63 

5.83 

1.52 

3.81 ? 

1.31 

1.98 ? 

5.98 

6.23 

5.51 

3.85 

4.17 

2.20 

2.68 

2.73 

1.35 



0.82 ? 
0.76 ? 
0.82 ? 



<W-fcalc> = °- 010 



<4alc//obs> = 0.865 



underestimated by a factor of 11 in our calculations and 
is not listed in table ITTI We notice that the other terms in 
the series 6d 5 D and 7d 5 D were not detected by SWBH. For 
example, the other multiplets in the Rydberg series, 3p P 2 - 
6d 5 D? j2 ,3 at A 4967.88 and 3p 5 P 2 -7d 5 D? i2j3 at A 4772.91 
should be theoretically stronger by 2.3 and 1.4 times the 
intensity of the A 4654.56 m ultiplet respectively a ccord ing to 
the recombination rates of lEscalante fc Victor! (jl992) , and 
should be above the noise in their respective parts of the 
spectrum. As noted in section 15.41 an identification as a 
Nil line is more likely, and therefore the line at A 4654.56 
is probably a misidentification as a O I line. 



5.9 Variation of the SED 

We compared results calculated with the original SED of 
MG and the new SED described in section 14.11 to gauge 



Lower-Upper 


A (A) 


-^rec / -^calc 


I 


-^obs 


3s 3 S°-4p 3 P § 


4368.2 


0.011 


12.30 


9.64 ? 


3p 3 P-7d 3 D° 


5275.1 


0.000 


6.34 


1.86 ? 


3p 3 P-8s 3 S° 


5298.9 


0.000 


5.01 


3.71 ? 


3p 5 P 3 -5d 5 D° 


5330.7 


1.000 


0.14 


0.23 ? 


3p 3 P-6d 3 D° 


5512.8 


0.001 


8.64 


3.13 ? 


3p 3 P-7s 3 S° 


5554.8 


0.000 


6.87 


1.59 ? 


2p 4 1 D 2 -2p 4 l So 


5577.3 


0.000 


1.64 


2.63 


3p 3 P-5d 3 D° 


5958.6 


0.003 


5.91 


4.75 ? 


3p 3 P-6s 3 S° 


6046.2 


0.000 


11.20 


8.81 ? 


3p 5 Pi-4d 5 D° 


6156.0 


1.000 


0.10 


0.18 ? 


3p 5 P 2 -4d 5 D° 


6156.7 


1.000 


0.17 


0.36 ? 


3p 5 P 3 -4d 5 D° 


6158.1 


1.000 


0.23 


0.42 ? 


2p 4 3 P 2 -2p 4 X D 2 


6300.3 


0.000 


146.00 


217.53 


2p 4 3 Pi-2p 4 X D 2 


6363.8 


0.000 


46.80 


75.94 


op r2 OS D 2 


6454 4 


1 ono 

l.UUU 


07 
u.u / 


n 1 n ? 


3p 5 P 3 -5s 5 Sg 


6456.0 


1.000 


0.10 


0.13 ? 


3p 3 P-4d 3 D° 


7001.9 


0.005 


10.80 


8.49 ? 


3p 3 P-5s 3 SJ 


7254.2 


0.002 


8.57 


15.64 ? 


3s 5 S°-3p 5 P 3 


7771.9 


1.000 


1.95 


3.52 


3s 5 S§-3p 5 P 2 


7774.2 


0.999 


1.39 


2.15 


3s 5 S°-3p 5 Pi 


7775.4 


1.000 


0.84 


1.30 


3s 3 SJ-3p 3 P 


8446.4 


0.017 


135.00 


114.19 


3p 5 Pi-3d 5 D° 


9260.8 


1.000 


0.31 


0.47 


3p 5 P 2 -3d 5 D° 


9262.6 


1.000 


0.52 


0.97 


3p 5 P 3 -3d 5 D° 


9265.8 


1.000 


0.73 


1.25 


{-^rec/-^calc) 


= 0.460 




(-^calc/-^obs) 


= 1.096 



§ Terms with no subscript are blended transitions of several 
states of angular momentum J. 



the dependence of the fluorescence excitation on theoretical 
assumptions about the stellar atmosphere. 

With few exceptions, most of the variation in results 
occurs in the Nil lines because of their high fluorescence 
excitation. The new SED increases the intensities of some 
Nil lines from p states by 25 per cent and the lines orig- 
inating in the 3d 3 D§ state by a factor of 2. The only 
significant decrease is a factor of 2 in the line 3p 3 Di- 
3d 3 Pg A 4459.937. We traced back the cause of those vari- 
ations to an increase of the stellar flux in the new SED at 
the wavelength of the transition 2p 2 3 P 2 -3d 3 D 2 A 533.815, 
and a decrease at the wavelength of the transition 2p 2 3 Pi- 
3d 3 PqA 529.413, which pump the observed lines. Minor vari- 
ations occur in Nil sing let lines Sp^i-Sd^ A 4447.03 and 
3p 1 D 2 -3d 1 F§A 6610.56, which are strongly excited by fluo- 
rescence. 



6 SUMMARY 

We have calculated the intensities of weak emission lines of 
Cn, Nl, Nil, Oi and On using a consistent nebular model 
and SED, which proved successful in previous work in repro- 
ducing the intense emission lines and overall density diag- 
nostics and surface brightness maps of IC 418. We compared 
our calculations with a published deep spectroscopic survey 
taking into account the fine structure of the lines and in- 
cluding some weak lines with possible blends. Although this 
nebula shows an almost spherical symmetry, our calculations 
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do take into account the position, size and orientation of the 
aperture used in the observations. 

As in our previous work on the Orion nebula, we found 
that fluorescence is responsible for the excitation of lines 
from s, p and d states in N II. Fluorescence also contributes 
significantly to the excitation of most lines of the s, p and d 
states in C n and some p and d states of O n while the rest 
of the 11 states are mostly excited by recombination. Some 
lines excited by fluorescence are sensitive diagnostics of the 
details of the SED, the expansion velocity gradient and tur- 
bulence of the nebula, the transition probabilities involved 
in the pumping mechanism, and their intensity is closely re- 
lated to lines in the stellar atmosphere. The recombination 
rates can explain most of the observed intensities from f and 
g states of all species in the ionized region although there are 
many observed transitions with excited core configurations 
awaiting the calculation of their effective recombination co- 
efficients. 

Ionic abundances deduced from dipole-allowed lines as- 
suming pure recombination rates without taking into ac- 
count fluorescence excitation of the lines can be overesti- 
mated by factors of up to 10 for Cll, 200 for Nil, 2 for 
O II, and 500 for N I and O I depending on the transition 
being considered, and a detailed nebular model and SED 
are needed to properly estimate ionic abundances from lines 
mostly excited by fluorescence. 

The automated procedure of line identification with the 
emili code seems to be in accordance with our calculations 
because the largest discrepancies between predicted and ob- 
served line intensities tend to occur in lines with dubious 
identifications or blends as indicated by that code. A deep 
spectroscopic nebular survey coupled with an appropriate 
line identification procedure can prove useful to test the ac- 
curacy of quantum mechanical calculations of atomic data. 

The N I quartet and O I triplet lines are produced by 
fluorescence in a narrow and dense interface of the neutral 
and ionized regions. Electron collisions compete with fluo- 
rescence in the excitation of the [N i] A 5200.26 and A 5197.90 
lines while fluorescence contributes a negligible amount to 
the excitation of the [O i] A 6300.30, 6363.78 and A 5577.34 
lines. However our calculations tend to underestimate sys- 
tematically the observed intensities of the [N i] lines by a 
factor of 5, and further work on the modeling of the PDR is 
necessary to fully understand the formation of those lines. 
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